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We study the extreme statistics of TV non-intersecting Brownian motions (vicious walkers) over a 
unit time interval in one dimension. Using path-integral techniques we compute exactly the joint 
distribution of the maximum M and of the time tm at which this maximum is reached. We focus in 
particular on non-intersecting Brownian bridges (" watermelons without wall" ) and non-intersecting 
Brownian excursions ("watermelons with a wall"). We discuss in detail the relationships between 
such vicious walkers models in watermelons configurations and stochastic growth models in curved 
geometry on the one hand and the directed polymer in a disordered medium (DPRM) with one 
free end-point on the other hand. We also check our results using numerical simulations of Dyson's 
Brownian motion and confront them with numerical simulations of the Polynuclear Growth Model 
(PNG) and of a model of DPRM on a discrete lattice. Some of the results presented here were 
announced in a recent letter [J. Rambeau and G. Schehr, Europhys. Lett. 91, 60006 (2010)]. 



I. INTRODUCTION AND MOTIVATIONS 

Extreme value statistics (EVS) is at the heart of op- 
timization problems and, as such, it plays a crucial role 
in the theory of complex and disordered systems [TJ [2] . 
For instance, to characterize the thermodynamical prop- 
erties of a disordered or glassy system at low temper- 
ature, one is often interested in computing its ground 
state, i.e. the configuration with lowest energy. Simi- 
larly, the low temperature dynamics of such systems is 
determined, at large times, by the largest energy barri- 
ers of the underlying free energy landscape. Thus the 
distribution of relaxation times is directly related to the 
statistics of the largest barriers in the systems. Hence, 
EVS has been widely studied in various models of statis- 
tical physics during the last few years [3]-[S] . This renewal 
of interest for EVS is certainly not restricted to physics 
but extends far beyond to biology H], finance [TUl |TT] 
or environmental sciences |12) where extreme events may 
have drastic consequences. 

EVS of a collection of M random variables 
Xi,--- ^Xj^f, which are identical and independent, or 
weakly correlated, is very well understood, thanks to 
the identification of three different universality classes in 
the large N (thermodynamical) limit [13]. In that case, 
the distribution of the maximum Xmax (or the minimum 
-''^min), properly shifted and scaled, converges to either 
(i) the Gumbel, or (ii) the Frechet or (iii) the WeibuU 
distribution, depending on the large argument behavior 
of the parent distribution of the XiS. By contrast, much 
less is known about the EVS of strongly correlated vari- 
ables. In this context, recent progress has been done 
in the study of EVS of one-dimensional stochastic pro- 
cesses which provide instances of sets of strongly corre- 
lated variables whose extreme statistics can be studied 
analytically. These include Brownian motion (BM) and 
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its variants ~ which can be often mapped onto elastic lines 
in 1 + 1 dimensions - [14H23j . continuous time random 
walks and Bessel processes [23] , or the random accelera- 
tion process (RAP) [211 [55]. Incidentally, it was shown 
that extreme value questions for such one-dimensional 
processes have nice applications in the study of random 
convex geometry in two dimensions. For instance, the 
extreme statistics of one-dimensional Brownian motion 
enters into the study of the convex hull of planar Brow- 
nian motion [22l [25 while extreme statistics of the RAP 
appears in the Sylvester's problem where one studies the 
probability p„ that n points randomly ch osen in the 
unit disc are the vertices of a convex n-sided polygon 
[27] . These stochastic processes involve a single degree 
of freedom, or several non-interacting degrees of freedom 
as in Ref . [201 1221 123] , where N independent Brownian 
motions are considered. A natural way to include the 
effects of interactions in a "minimal", albeit non trivial, 
model is to constraint these N Brownian motions not to 
cross: this yields a model of N non-intersecting Brownian 
motions, which we will focus on in this paper. 

Here we thus consider TV non-colliding Brownian mo- 
tions 

xi(r) < ••• <a:jv(T) , Vre [0,1] , (1) 

on the unit time interval, r e [0,1]. We will consider 
different types of configurations of such vicious walkers, 
which are conveniently represented in the {x, t) plane. 
In the first case, called "watermelons", all the walkers 
start, at time t — Q and end, at time t = 1, at the 
origin (see Fig. |4]). They thus correspond to TV non- 
intersecting Brownian bridges. In the second case, we 
will consider such watermelons with the additional con- 
straint that the positions of the walkers have to stay 
positive (see Fig. [?]): we will call these configurations 
"watermelons with a wall" and they thus correspond to 
non-intersecting Brownian excursions. Finally, we will 
also consider "stars" configuration, where the endpoints 
of the Brownian motions in r = 1 are free. They corre- 
spond to non-intersecting free Brownian motions, a real- 
ization of which is shown in Fig. |6]for TV = 2. 



Motivated by extreme value questions, we will compute 
here the joint distribution Pisf{M,TM) of the maximal 
height and the time at which this maximum is reached 
(see Fig.|4]): 

M = max Xn{t) , xn{tm) = M . (2) 

0<T<1 

The (marginal) distribution of M has recently been stud- 
ied by several authors pMM| . while we have announced 
exact results for this joint distribution in the two first 
cases (for bridges and excursions) in a recent Letter [35] . 
The goal of this paper is to give a detailed account of the 
method and the computations leading to these results. 
On this route we will also provide several new results, 
including for instance results for the star configuration. 

In the physics literature, such non-intersecting Brow- 
nian motions were first introduced by de Gennes in the 
context of fibrous polymers |36j . They were then widely 
studied after the seminal work of Fisher [37] , who named 
them "vicious walkers", as the process is killed if two 
paths cross each other. These models have indeed found 
many applications in statistical physics, ranging from 
wetting and melting transitions [37], commensurate- 
incommensurate transitions |38j . networks of poly- 
mers [55] to persistence properties in nonequilibrium sys- 
tems [4cr. 

Vicious walkers have also very interesting connections 
with random matrix theory (RMT), in particular through 
Dyson's Brownian motion [?T]. For example, for water- 
melons configurations, it can be shown that the positions 
of the N random walkers at a fixed time r, correctly 
scaled by a t dependent factor, are distributed like the 
N eigenvalues of the random N x N Hermitian matri- 
ces of the Gaussian Unitary Ensemble of RMT corre- 
sponding to /3 = 2 [32] . If one denotes by Pjoint (x, t) = 
-Pjoint(2;i, • • • ,xn,t) the joint distribution of the posi- 
tions of the walkers at a given time r, in the watermelons 
configuration, one has indeed (see for instance Ref. [30] 
for a rather straightforward derivation of this result): 

Pjoi„t(x,r) = Z^V(r)-~' Jl {x, - x,fe~^) ,(3) 

l<i<j<N 

where Zj^ is a normalization constant and (j{t) — 
r(l — r) and where we use the notation x'^ — 
Tlif^i^l- This means in particular for the top path 
that xn{t) / \/t{\ — t) is, at fixed r, distributed like the 
largest eigenvalue of GUE random matrices, which means 
that in the large N limit one has [33] 

where X2 is distributed according to the Tracy- Widom 
(TW) distribution for (i — 2, F2, namely Pr[x2 < x] — 
T2{x). In the case of watermelons with a wall, the po- 
sitions of the random walkers are instead related to the 
eigenvalues of Wishart matrices [301 133] and the fluctua- 
tions of the top path x n (r) are again described by J-2 (x) . 
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FIG. 1: Dynamical rules of the PNG model (in continuous 
time) as described in the text. 



Yet another reason why there is currently a rekin- 
dled interest in vicious walkers problems is because of 
their connection with stochastic growth processes in the 
Kardar-Parisi-Zhang (KPZ) universality class [45l [46] in 
1 + 1 dimensions. This connection is believed to hold 
for any systems belonging to the KPZ universality class 
|48) . and it can be rigorously shown on one particu- 
lar model of stochastic growth, the so-called polynuclear 
growth model (PNG) [31 [SO]- It is defined as follows 
(see Fig. [T]). At time t = a single island starts spread- 
ing on a flat substrate at the origin a; = with unit 
velocity. Seeds of negligible size then nucleate randomly 
at a constant rate p = 2 per unit length and unit time 
and then grow laterally also at unit velocity. When two 
islands on the same layer meet they coalesce. Mean- 
while, nucleations continuously generate additional lay- 
ers. In the flat geometry nucleations can occur at any 
point X while in the droplet geometry, which we will fo- 
cus on, nucleations only occur above previously formed 
layers. Therefore, for the droplet geometry, denoting by 
^drop(a;j t) the height of the interface at point x and time 
t, one has /idrop(a;,i) = for |a;| > t. On the other 
hand, in the long time limit, the profile for \x\ < t 
becomes droplet-like (/idrop(a;, i)) ~ 2t-\/l — {x/tY [51j . 
but there remain height fluctuations around this mean 
value (see Fig. [2] left). The standard way to characterize 
these fluctuations is to look at the width of the interface 
WL{t) = ((/idrop(a;,i) - {K,op{x,t))f)h ^ L'^Wit/L-) 
with universal exponents C = 5 and z = | [52| I53j . 
in agreement with the fact that this model belongs to 
the KPZ universality class. More recently, it was shown 
for several models belonging to the KPZ universality 
class [54H59] that, in the growth regime to <^ t <^ 
(where tg is a microscopic time scale) universality extends 
far beyond the exponents ^ and z but also applies to full 
distribution functions of physical observables. In partic- 
ular, the scaled cumulative distribution of the height field 
at a given point coincides with the TW distribution J"^ 
with /3 = 2 (respectively /3 = 1) for the curved geom- 
etry (respectively for the flat one), which describes the 
edge of the spectrum of random matrices in the Gaussian 
Unitary Ensemble (respectively of the Gaussian Orthog- 
onal Ensemble) [33] • Height fluctuations were measured 
in experiments, both in planar |62j and more recently in 
curved geometry in the electroconvection of nematic liq- 
uid crystals [63] and a good quantitative agreement with 
TW distributions was found. 
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Although the relation with random matrices were ini- 
tially achieved |54j through the longest increasing sub- 
sequence of random permutations [SHI HI] , using in par- 
ticular the results of the seminal paper by Baik, Deift 
and Johansson [60j, it was then realized that the PNG 
model in the droplet geometry is actually directly related 
to the vicious walkers problem in the watermelon config- 
uration |55| . This was shown through an extension of 
the PNG model to the so called multi-layer PNG model, 
where non-intersecting paths naturally appear. Indeed, 
one can show that the fluctuations of the height field of 
the PNG model in the droplet geometry hdiop{x,t) are 
related, in the large time limit to the fluctuations of the 
top path XNir) for the vicious walkers problem in the wa- 
termelon geometry in the large N limit. This mapping, 
for N,t:$> 1 reads [SSIIM] : 




~A2{u)~u^, (5) 

where A2{u) is the Airy2 process [55,^ which is a sta- 
tionary, and non-Markovian, process. In particular, 
Pr[^2(0) < x] — J^2{x), which is consistent with Eq. Q. 
Hence, from Eq. ^ xn and r map onto h and x in 
the growth model while N plays essentially the role of 
t (to make the correspondence between N and t exact 
one has to consider a watermelon configuration in the 



interval t G [0, iV]). To characterize the fluctuations of 
the height profile in the droplet geometry beyond the 
standard roughness WL{t) it is natural to consider the 
maximal height M and its position Xm [3S] (see Fig.l2|. 
According to KPZ scaling, one expects M — 2t ^ r^'^ 
while Xm ^ t^^'^- On the other hand from Eq. ([H]), 
the joint distribution Pt(M, of M and Xf,[ can be 
written as 

Pt{M,XM)^t-^VAuy({M-2t)r-3,XMr"^) , (6) 

where 7^Airy(2/j a;) is the joint distribution of the maxi- 
mum y and its position x for the process A2(u) — u^. 
Finally, the relation ^ also gives us some interesting in- 
formation for the vicious walkers problem. In the large 
N limit, one has indeed from Eq. ([5| that the joint dis- 
tribution Piq{M,TM) for the watermelons configuration 
is also given by 'PAiry (y, a^)- One has indeed, for 1: 

Pn{M,tm) 




and therefore Eq. ^ and Eq. ([t]) show that one can 
obtain the distribution of Pt{M,XM) for the growth 
model from the large N limit of the joint distribution 
P]\j(M,tm), which we compute here for any finite N. 



h{x,t) 




FIG. 2: Left : Height profile hdiop{x,t) at fixed time t, as a function of x for the PNG model in the droplet geometry. Xm is 
the position at which the maximal height M is reached. Right : Directed polymer with one free end. Here M corresponds 
to the energy of the optimal polymer, and M — 2t ~ 0{t^^'^) while Xm ~ 0{t^^'^) corresponds to the transverse coordinate of 
the end point of this optimal polymer. 

I 



It is also well known that stochastic growth models in [351 137] . This is also the case of the PNG model and to 
the KPZ universality class can be mapped onto the model make this connection as clear as possible, we consider 
of the directed polymer in a disordered medium (DPRM) a discrete version of the PNG model introduced by Jo- 
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FIG. 3: Left ; Model of a directed polymer in a random medium corresponding to the PNG in the droplet geometry ([lOf . On 
each site • there is a random energy variable w{i,j) which corresponds to a nucleation in the growth model. In red and blue 
we have drawn two configurations of the polymer which belong to the ensemble ^Pt(0, X). Right : Model of a directed polymer 
in a random medium corresponding to the PNG in the flat geometry (TTJ. In red and blue we have drawn two configurations 
of the polymer which belong to the ensemble 0) which is equivalent to ^t(0, X) shown on the left panel. 



hansson [BS]. It is a growth model with discrete space 
and discrete time. The height function h{x, t) is now an 
integer value h{x,t) e N while a; S Z and t e N. The 
dynamics is defined as follows 

/i(a;,t + l) = m.ayi[h[x — l,t)^h{x,t)^h{x + 1,1)] 

+ w(x,t + l) , (8) 



where the first term reproduces the lateral expansion of 
the islands (and also their coalescence when two of them 
meet) in the continuous PNG model while ijj{x, <) S N is a 
random variable, distributed independently from site to 
site, which corresponds to the random nucleations (see 
Fig. [1]). 



In the droplet geometry, one has ijj{x,t) = if |a;| > t 
while this constraint is removed in the flat geometry. To 
understand better the connection between this model Q 
and a model of directed polymer, we perform a rotation 
of the axis {x,t) (Fig. Isl) and define 



^{i- i - 1) = w{i,j) 



(9) 



One can then simply check (for instance by induction on 
t), that the height field /idrop(2^,0 evolving with Eq. ([s]) 
in the droplet geometry, is given by 



hdtoniXA) = max 

ce<Pt(o,x) 



(10) 



where *Pf (0, X) is the ensemble of directed paths of length 
i, starting in a; = and terminating in X (Fig. |3]). 



This formula ( 10 ) establishes a direct link between the 
PNG model and a model of directed polymer. Indeed, 
the height field /idrop [X, t) corresponds the energy of the 
optimal polymer (which is here the polymer with the 
highest energy) with both fixed ends in and X and 
where the random energies, on each site, are given by 
w{i,j). The aforementioned result for the distribution of 
hdrop{X,t) shows that the fluctuations of the energy of 



the optimal polymer are described in this case by J-2, the 
Tracy- Widom distribution for /? = 2. The study of this 
distribution, for the DPRM in continuum space, has re- 
cently been the subject of several works, both in physics 
|66 | and in mathematics j67j . 

Similarly, one can also write the height field /iflat(0,i) 
in the flat geometry (here of course the fluctuations of 
the height fleld are invariant by a translation along the 
X axis). One has indeed : 



/iflat(0,i)= max 

-t<x<t 



max 



(»j)ec 



, (11) 



where *Pt {X, 0) is the ensemble of directed paths of length 
t, starting in X and ending in (see Fig.lsl). Given that 
these two ensembles ^t{X, 0) and ^t(0, are obviously 
similar (Fig. [s]), these two equations (10 1 and (11) allow 
to write [551 



Iflat 



(0,i) 



h. 



-t<x<t 



drop 



{X,t). 



(12) 



For the directed polymer, /iflat(0,i) corresponds to the 
energy of the optimal polymer, of length t and with one 
free end. From that Eq. (12 1, one sees that the joint dis- 



tribution of M and tm in the vicious walkers problem 
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corresponds to the joint distribution of the energy M 
and the transverse coordinate Xm of the free end of the 
optimal polymer. From the aforementioned results for 
the distribution of /iflat(0,i), we also conclude that the 
cumulative distribution of M is given by J^i, the Tracy- 
Widom distribution for /? = 1. Although this result was 
obtained rather indirectly using this relation (12), it was 



recently shown directly by computing the cumulative dis- 
tribution of the maximal height of N non-intersecting ex- 
cursions in the limit ^ oo, in Ref. [69] . We conclude 
this paragraph about the DPRM by noticing that this 
model, with one free end, has been widely studi ed in 
the context of disordered elastic systems [7D] , in particu- 
lar using the approximation of the so-called "toy- model" 
where the Airy2 process A2{u) in ([s]) is replaced by a 
Brownian motion [51 [TO] . A recent discussion of the toy 
model and its applications to the DPRM can be found 
in Ref. [71]. We also mention that physical observables 
related to Xm have been studied for the DPRM in dif- 
ferent geometry, like the winding number of the optimal 
polymer on a cylinder |72[ . 

The paper is organized as follows. In section II we 
present in detail the method to compute the joint distri- 
bution of M and tm for N = 2 non-intersecting Brownian 
motions. We will show results for bridges, excursions and 
stars. In section III, we will extend this approach for any 
number N of vicious walkers both for bridges and excur- 
sions. We will then confront, in section IV, our analytical 
results to numerical simulations, which were obtained by 
simulating Dyson's Brownian motion, before we conclude 
in section V. Some technical details have been left in the 
appendices. 



(ii) free boundary conditions where we consider non- 
intersecting free Brownian motions, also called 
"stars" configurations (see Fig.|6|: this will be treated 
in subsection F, 

(iii) periodic boundary conditions, plus a positivity con- 
straint where we consider non-intersecting excursions, 
also called "watermelons with a wall" configurations 
(see Fig. [S]): this will be treated in subsection G. 

We compute, for each case, the joint distribution of the 
couple (M, Tjv/) which are respectively the maximum M 
of the uppermost vicious Brownian path, and its position 
Tm- M = maxo<r<i a;2(T) = 2^2 (tm)- Integrating over 
M, we find the probability distribution function (pdf) of 
the time to reach the maximum tm- 

This computation is based on a path integral approach, 
using the link between the problem of vicious Brownian 
paths and the quantum mechanics of fermions. This is 
described in subsection B. Then we present, in subsection 
C, the method to compute the joint probability distribu- 
tion function of the maximum and the time to reach it 
for vicious Brownian paths, with general boundary condi- 
tions. An emphasis is put on the regularization scheme, 
essential to circumvent the basic problem of Brownian 
motion which is continuous both in space and time: you 
cannot force it to be in a point without re-visiting it in- 
finitely many times immediately after. At this stage, gen- 
eral formulae will be obtained, and will be used directly 
in the paragraphs E, F and G. The extension to a general 
N number of paths will be treated straightforwardly in 
Section III. 



II. DISTRIBUTION OF THE POSITION OF 
THE MAXIMUM FOR TWO VICIOUS 
BROWNIAN PATHS 

A. Introduction 

Our purpose is to study N simultaneous Brownian mo- 
tions Xi{t), X2{t) . . . ,xi^{t), subjected to the condition 
that they do not cross each other. In the literature, such 
system has been studied first in the continuous case by de 
Gennes [35], and in the discrete case by Fisher [37] ■ He 
named this kind of system vicious walkers, as the process 
is killed if two paths cross each others. The name 'vicious 
Brownian paths' will refer in the following to Brownian 
paths under the condition that they do not cross each 
other. 

In this Section, we focus on the case N — 2. The 
two vicious Brownian paths start at 2:1(0) — 2:2(0) — 
and obey the non-crossing condition xi(t) < 2:2 (t) for 
< T < 1. We examine three different geometries: 

(i) periodic boundary conditions where we consider non- 
intersecting bridges, also called "watermelons" con- 
figurations (see Fig. |4]) : this will be treated in sub- 
section E, 



B. Path integral approach: treating the vicious 
Brownian paths as fermions 

Let us consider two Brownian motions, obeying the 
Langevin equations 



(J = l,2), 



(13) 



where ^j(r)'s are independent and identical Gaussian 
white noises of zero mean, i.e. (Cj (''')) = and 
(Cj(''')Ci ('''')) = Taking into account the non- 

crossing condition, the probability that the two paths end 
at (xi(tb) = 61,2:2(^6) = 62) at time tf,, given that they 
started at {xi(ta) = 01,2:2(^0) — 0,2) at time ta is 

P2{bi,b2,tb\ai,a2,ta) = ^ Jvxiir) J'Dx2{t) 



tb l/d^lM 
2 V d 



Xl{ta} = 0'i X2{ta) = a2 

n 0[x2iT)-x,ir)]], (14) 

ta<T<tb 
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where the subscript '2' in the notation p2 means that 
we treat the case of iV = 2 Brownian paths. The ex- 
ponential terms are the weight of each Brownian path, 
and the Heaviside theta function ensures that the paths 
do not cross. In this Section, bold letters will denote 
vectors with two components, such as for example the 
positions at a given time of the two Brownian paths 
x{t) = (xi(t),X2(t)) or the starting points a = (01,02). 
In the previous formula ( 14 ) , z is a normalization con- 
stant, such that 



+00 



d&2 



62 



d6i p2{b,tb\a,ta) 



(15) 



Notice that the integration is ordered, and we will write 
it shortly /„^,db^/+rd62/!:,d6i. 

We recognize in formula ( [l4| the path integral repre- 
sentation of the quantum propagator of two identical free 
particles, subjected not to cross each other. In one di- 
mension, this can be implemented by requiring that the 
two particles are fermions |30l [36] . Therefore the prob- 
ability P2{b, ti,\a, to) reduces to the quantum propagator 
of two identical fermions in one dimension: 



P2[b,h\a,ta) = {b\e 



(16) 



with iJfrcc the total Hamiltonian of the two-particle sys- 
tem TJfroc = i?/,eo ® I^^^ + I^^^ ® ^frol' being the iden- 
tity operator acting in the Hilbert space of the j*^ par- 
ticle. The one-particle Hamiltonians are ^ff/oc — "IdiT^ 
(for j = 1,2) in a position representation. The states 
\ai,a2) (and 161,62)) are the tensorial product of the 
one-particle eigenstates of position |ai,a2) = |ai) (g) [02), 
where each state \aj) obeys the eigenvalue equation 
\aj) = a-i\aj), with X^^^ the position operator act- 
ing in the Hilbert space of the j**^ particle, and Uj the 
position of this particle. 

In our calculations we deal with more general Hamil- 
tonians, but they will be written in the same way as 
the sum of individual one-particle Hamiltonians H = 
(g) + 1(1) (g, ij(2) (no interaction between the two 
particles), and with H'^^\x,p) — H^'^\x,p) (identical 
particles). Thus we will treat systems of two identical 
indistinguishable particles, so that one can write down a 
two-particle eigenstate of the total Hamiltonian H as an 
antisymmetric product of the one-particle eigenstates of 
H^^\ Explicitely, if \Ei) and |i?2) are two eigenstates of 
the one-particle Hamiltonian H'^^\ with energies Ei and 
E2, then the eigenstate of the two-particle problem with 
energy E = Ei + E2 is 

4f (|Si) ® \E2) - \E2) ® \E{)) . (17) 



\E) = ^ (|i?i) ® \E2) - \E2) ® \E,)) 

Projecting onto the position basis \x) = \xi) ® \x2), and 
using the general notation (pE^ixj) = {xi\Ej) (lower-case 
letter (}> for one particle wave function), and ^e{x) = 



{x\E) (capital letter $ for the two particle wave func- 
tion), this is the usual Slater determinant 



'^e{xi,X2) 



2! 



5i {xi)4>E2 {X2) - (t>E2 {xi)(t>Ex {X2)) 



-;= det d)E(xi). 
/2I i<i.j<2 ■" 



(18) 



Hence inserting the closure relation I = J2e\-E)(-E\ 
into (16), one obtains the spectral decomposition of the 



propagator (with E ^ Ei + E2) 

P2{b,tt\a,ta) = J2 <i>Em*Eia)e-^'''+^^^^'^~'-\ 

El ,£'2 

(19) 

This result can also be found via the Karlin-McGregor 
formula [75] . It states that the probability of propagating 
without crossings is nothing but the determinant formed 
from one-particle propagators (named pi ) : 

P2{b,tt\a,ta) = det pi{bi,tt\aj,ta). (20) 

l<ij<2 



The discrete version of this formula ( 20 ) is given by the 



Lindstrom- Gessel-Viennot (LGV) theorem [71]. The ba- 
sic idea behind this formula is nothing but the method of 
images. To recover our quantum mechanical expression 
obtained by considering fermions ( 19 ), one has to express 
the one-particle propagator as 



pi{bi,tt\aj,ta) ^J2^Eib,)4>*E{aj)e 



-ta) 



(21) 



and use the Cauchy-Binet identity (B7| 



C. Joint probability distribution of position and 
time 

The method to compute the distribution of the maxi- 
mum and the time at which it is reached can be described 
as follows. Let us consider the process {xi{t),X2{t)) 
of such Brownian paths with the non-intersecting con- 
dition with ta = and tf, — T. Keeping the same nota- 
tions, the starting points are a;i(0) = ai and a;2(0) = 02, 
or a;(0) = a, and the end-points are Xi{T) = 61 and 
X2{T) ~ 62, or x{T) — b. First we consider a and b 
as fixed. Given a set of points y — (2/1,2/2), and an 
intermediate time Tcut, such as < Tcut < T, we ask 
what is the probability that 2;i(Tcut) G [yiiUi + dj/i] 
and X2iTcut) e [2/2,2/2 + d2/2] ? In other terms, we 
seek Q2{y,Tcut\b,a^T) the joint probability density of 
(y,rcut), given that the paths start in a at r = and 
end in 6 at r = T. The answer to this question is well 
known from quantum mechanics: cut the process in two 
time-intervals, one from t = to r = Tcut, and the second 
piece from r = Tcut to r = T. Because of the Markov 
property, these two parts are statistically independent, 
and one should take the product of propagators to ob- 
tain the joint pdf of (y,Tcut) = (2/1, 2/2, Tcut), given that 
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the initial and final conditions are respectively (a, 0), and 
(6,T): 



P2ib,T\y,Tcnt) X P2(j/,Tcut|a,0) 

^2(b,a,T) 



<92(y,Tcut|b,a,T) = 

(22) 

where Z2{b, a, T) is the normalization constant, obtained 
by integrating over all possible cuts Tcut , and all possible 
intermediate points 



Z2= j drcut / dy P2(b, T|y, Tcut) x P2(y, Tcut|a, 0) , 
J J orcl 



(23) 

where we remind that the integration over y is ordered: 
dy = dj/2 d?/i . Using the closure relation of 
the states |y), one finds that 



Z2(b,a,r) =Tp2(b,r|a,0) 



(24) 



Because we will not care about the lowest path xi(t) 
and its intermediate position j/i in the following, we com- 
pute the marginal by integrating over all possible values 
of yi = xi(rcut), obtaining P2(2/2, Tcutl^, a, T) the joint 
probability density that the upper path is in ?/2 at time 
Tcut given that the paths start in a at time t = and 
end in b at time t = T: 

dyi Q2(y,Tcut|b,a,T) 

oo 

j-^'^Vl P2(b,T|y,7-cut) X P2(y,Tcut|a,0) 



Tp2{b,T\a,Q) 



(25) 



From now on we will set, for simplicity, T — 1. This can 
be done without loss of generality because the Brownian 
scaling implies M cx T^/^ and tm oc T. A simple di- 
mensional analysis allows to resinsert the time T in our 
calculations. 



D. Maximum, and the regularization procedure 

The maximum M and the time to reach it tm are 
defined by 



M = max X2{t) = X2{tm) 

0<T<1 



(26) 



Hence, to compute the joint probability of the couple 
(M, Tjvf), we want to impose y2 = M and Tcut = tm^ 
knowing that the upper path does not cross the line 
X = M. This condition is implemented by inserting the 
product of Heaviside step functions 



n 6[M-X2{t)], 



(27) 



0<r<l 



in the path integral formula (14 1. The propagator of two 



vicious Brownian paths, with the constraint that they do 



not cross x = reads 



P<AL2ib,'tb\a,ta) = 



Z<M 



Vxi{t) [ Vx2{t) 



a;i(ta)=ai 



{ [] e{M^X2{T))6{M-X^{T))e{x2{T)-Xl{T)) 



ta<r<tt 



X e 



(28) 



for all < ta < < 1, ai < a2 < M and 6i < 62 < M. 
The second Heaviside step function does not change the 
value of the path integral (it does not select a smaller 
class of paths) because for all ta < t < tb, if xi{t) < 
X2{t) and 2:2 (r) < M, then xi{t) < M automatically. 
Accordingly to our first analysis, this propagator (28) 



reduces to the propagator of two identical fermions with- 
out interaction between them, but with a hard wall in 
X = M . Hence the associated Hamiltonian is 



H<iM — H_ 



(1) 



^(2) 



(2) 
<M 



with the one particle Hamiltonian 

1 d2 



H 



<M 



2dx 



2 + V<M{xj), 



(29a) 



(29b) 



for j — 1,2, where V^m{x) is the hard wall potential in 
X = M: 



V^m{x) 







if -00 <x < M 



-00 if X > M . 
This one particle Hamiltonian has eigenfunctions 



(l^kAxj) 



with the associated energies 



Ek, - kj > 



(29c) 



(29d) 



(29e) 



Hence, with these informations, one is able to compute 
the propagator p^M{b, ti,\a,ta) using the spectral decom- 
position, as in Eq. (19): 



P<M,2i^,ti,\a,ta) = (b|e 



-{th — ia)H<M I 



dfci / dk2^k{bWk{a)e 



-(tb-ta)^- 



(29f) 



with $fc(a;) the Slater determinant built with the one 
particle eigenfunctions (29d). 



To obtain the joint pdf of M and ta/, one would insert 
this propagator in formulas ( 22p5 ), with y2 — M and 
Tcut = Tm- With this procedure one naively finds zero, 
because we impose the intermediate point 1/2 to be on 
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the edge of the hard wall potential. This problem origi- 
nates in the continuous nature of the Brownian motion: 
once in a point, the Brownian motion explore its vicinity 
immediatly before and after |75j . This implies that we 
cannot prescribe the upper path to be in X2{tm) = M 
without being in X2(t) > M for t close to tm- Another 
problem is that the normalization constant can not be 
written as a single propagator. Indeed one can not insert 
the closure relation as we did from Eq. ( 23 ) to Eq. ( 24 1 , 



because the potential depends actually on the position 
2/2 = M. 

Hence one should prescribe a regularization scheme to 
avoid this phenomenon. A common way to deal with it 
is to take j/2 ~ 2^2 (''"Af ) = M — rj, with 77 > a small regu- 
larization parameter, compute all quantities as functions 
of T], and at the end take the limit 77 — >■ 0. This regular- 
ization scheme is similar to that introduced in Ref. [TO] . 
Following the same ideas, one obtains the joint pdf as 

P2 (M, Tm b, a = hm , 30a) 

»?^o Z2[r]\b,a) 

where W2{M — ri,TM\b,a) is the probability weight of 
all paths starting in a and ending in b in the unit time 
interval and such that X2{tm) = M — rj, and Z2{r]\b,a) 
is the normalization constant. Explicitly, one has 



y2 



W2{M-r,,TM\b,a)^ J dy2 J dyi {s{y2~ [M ~r])) 

— 00 —00 

X P<M,2ib,l\y,TM) P<Af,2(y,TM|a,0)|, (30b) 

for all M > max{a2,62}, and where we write a dummy 
integration over y2 with the delta function forcing 7/2 to 
be M — r). The normalization depends on the regulariza- 
tion parameter 77 

1 °° 

Z2{T]\b,a)= j dTM j dM W2{M -r^,TM\h,a) . 

max{a2 ,62 } 

(30c) 

In the following paragraphs, all three configurations 
have their paths beginning at the origin. The problem of 
the continuous Brownian paths forbids to take directly 
a = (0,0). Instead of that, we will separate artificially 
the paths by an amount e > 0. Details of this regulariza- 
tion are left in the concerned paragraphs. 



E. Periodic boundary condition: N = 2 vicious 
Brownian bridges 

In this case, the two paths start from the origin at 
r = and arrive also at the origin at final time t — 1: 
xi(0) = 0:2(0) = a;i(l) = a;2(l) = 0. Without extra 
condition, this defines two Brownian bridges. Here we 
compute the joint probability density function for two 



x(x)' 
M 










Vv. 





FIG. 4: One realization of A'^ = 2 non-intersecting bridges, 
i.e. a watermelon configuration, with the maximum M and 
the time tm at which this maximum is reached. 



Brownian bridges, under the condition that they do not 
cross each other. 

As discussed before, there is again a regularization 
scheme to adopt in both starting and ending points. We 
separate by an amount e the two paths at the starting 
and ending points, taking a = 6 = (0, e) = e. At the 
end of the computation, we shall take the limit e — > 0. 
Then the joint probability density function for two vi- 
cious Brownian bridges is 



P2,B (M, TM ) = hm P2 (M, TM I e, e) 

e— >-0 



(31a) 



where P2 {M, t mI^, b ) is defined in Eq. ( 30a I . The nu- 
merator of Eq. (30a), the probability weight, reads 



W2,b{M ~Tj,TM,e) ^W2{M -r,,TM\e,e) , (31b) 



and the denominator of Eq. (30a), i.e. the normalization, 
is given by 



Z2,b{v,<^) = Z2{il\e,e), 



(31c) 



where the subscript 'B' refers to 'Bridges', and with the 
relation deduced from Eq. (30a) 



P2,b[M, Tm) = hm hm — r . (31d) 

We obtain this limit by expanding the numerator and the 
denominator in power series of 77 and e, keeping only the 
dominant term in each case. 

To compute the proba bility weight W2,b{M — rj, tm, e), 
we start with Eq. (30b), in which we put a = b = e: 



W2.B{M-r^,TM,e)^ / dy {s{y2-{M-rj)) 

J Ol d 

X p<A/,2(e, l|y,TM) P<A./,2(y,TA/|e, 0)| . (32) 
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Using the spectral decomposition ( 29f ) , one has 



W^2,B(A/-ry,TM,e) - / dy f %2 - (M-77)) 



X dfc $fc(y)$*(e)e-^-V| . (33) 

The integral over 2/2 is easy due to the delta function. By 
antisymmetry of Slater determinants in the exchange of 

k\ and k2, the product $fc(y)$fc(e)e~'^*'^ ~ is symmetric 
in the exchange of ki and k2- Then one can replace ^k{y) 
by twice the product of its diagonal terms which is 

^^ ^.,(,,)0..(A/-.) ^ A,^(^^)^^(^3). (34) 
v2' V 

The same operations applies to ^f,'{y). 

.4 rM-n 



X / dfc' k'^ $fc'(e)e 
/o 



dfc fc2 0fei(2/i) *fc(e)e"""''^ ■ (35) 



Permuting the order of integrations, one can compute the 
integration with respect to yi as 



M-rj 



dyi {yi)^k, iyi) = Sik, - k[) + 0{r^) , (36) 



due to the orthonormalization of the eigenfunctions 
(j)k{x). Performing the integration over k'l leads to, at 
lowest order in 77, 



00 00 



W2,BiM -Tj,TM,e) - J dfcie-T J d/c2 J dk'^ 



X <^ /c^,fc2$fei,fe'(e)$fc,,fc,(e)e 



+ 0{v^). (37) 

The next step is to expand in powers of e inside the Slater 
determinants: the first column does not depend on e, 
and in the second column, we expand at order e, each of 
the two elements, sweeping the constant term by linear 
combination with the first column. This yields 



'?!'fc2(0) '/'fe2(e) 
j_ 0fc,(O) (/-^^(O) 

=£—^62(91, 92) + 0(e') (38) 



where we use the scaled variables qi = kiM and the de- 
terminant 



©2(91, 92) = , det I g: 



cos [Qi - j- 



(39) 



Performing the same expansion in powers of e in 
'J'fei.fe^(e): with q'2 = k'2M, one can identify the domi- 
nant term as a factor of the product of powers rj'^e^, so 
that 

W^2.B(M-r/,TM,e) =r]VW2,s(Af,TAf)-l-0(7yV,r;2e3) , 

(40) 



Inserting the expansion (38) in (37), and identifying the 



the leading term, one obtains 

W2,i3(M,TM) 



dg2 92 ©2 (91,92)6 



X / dg2 9262(91,92)6 
/o 



(41) 



It can be written in the more compact form 
W2,b{M,tm) 



T2(gi|Af,l - TM)T2{qi\M,TM), (42) 
with the help of the function T2 



T2 (9i|A/,tm) 



dg2 926 62(91, 92) , (43) 



which can also be written as a determinant: 



T2 (9i|M,tm) = V 
cos (91 



M 

2 ) ^91 COS (91 - tt) 

M 



(44) 



Hi{x) and H2{x) are the first and second Hcrmitc poly- 
nomials. This expression is obtained by factorizing all 
the 92-dependence in the last line of the determinant, 
and performing the integration directly in this last line. 
The Hermite polynomials then appear naturally, see Ap- 
pendix [Ej Then one proceeds to the expansion of the de- 
terminants T2 with respect to their last lines, which per- 
mits the integration over 91 (four terms, because there are 
two determinants 2x2), and one obtains a simple, though 
long, expression (apart from the exponentials, only alge- 
braic terms appear). 



The limit in Eq. (31d) exists provided the normaliza- 



tion constant admits the following expansion 

Z2A'n, 6) = rfe^K2,B + oi-n'^) , (45) 
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FIG. 5: Distribution of the time at which the maximum is 
reached for N = 2 non-intersecting bridges. The solid line 
corresponds to our analytical formula in Eq. ( 48 1 while the 



symbols correspond to the results of our numerical simula- 
tions. 



where K2,b is a number (independent of M and tm)- It 
can be computed using the normalization of the joint pdf 



1 = 



/ dM / dTM P2,b{M,tm) 
Jo Jo 



. dM / dTM W2.BiM,TM) 
J^2.B Jo Jo 

One finds K2^b — 1/tt. The joint pdf reads 



(46) 



P2,b{M,tm) 



1 



^[TAf(l-TMF/2 

M 



Ho 




l-le-''''H2{V2M)] I (47) 



where for compactness we use the Hermite polynomial 
H2{x) = Ax^ — 2. A direct integration over M gives the 
probability distribution function of tm, the time at which 
the maximum is reached, regardless to the value of the 
maximum: 



P2.b{tm) = 41 



l + lOTMjl-TM) 
[l + iTMil-TM] 



15/2 



(48) 



The fact that tm enters in this expression only through 
the product tm(1 — tm) reflects the symmetry of the 
distribution around tj\/ = 1/2, which of course is ex- 
pected for periodic boundary conditions. To center the 



distribution, over an unit length interval, one can use 
Tm = 1/2 + um/'2, and the distribution of is then 



Pr 



centered / \ o 
2,B [^M) = 2 



(49) 



(2-4,)3/2 (2-4,)5/: 
Two asymptotic analysis can be made: 

• for Tm — (or equivalently tm — 1), one has 

P2,b{tm) - 120tIj + 0{tIj), 

• and for r^v/ in the vicinity of 1/2, better written in 
terms of the behaviour in um = of the centered 
distribution 



pccntorcd/ \ 
P2.B K^M) 



Furthermore, an integration of P2,b{M,tm) with re- 
spect to Tm gives the pdf of the maximum g{M) 
(the derivative of the cumulative distribution). Using 
the change of variables tm = (1 + sin((/3))/2, the integral 
gives 

FiBiM) = I dTM P2.b(M,tm) 
Jo 

which coincides with the result computed in Ref. [30] . 

F2^Bix) ^l-4x'^e-' ' 



F. Free boundary conditions: N = 2 stars 
configuration 



x(%)' 
M 








4 









WW 




/X 



FIG. 6: One realization of A'^ = 2 non-intersecting free Brow- 
nian motions, i.e. in "star" configuration, with the maximum 
M and the time tm at which this maximum is reached. 

In this subsection, we consider two vicious Brownian 
paths starting from the origin with free end-points: this is 
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the star configuration (a subscript 'S" will be used). The 
starting points are the same as before a = e = (0, e), but 
the endpoints are free, so that for M > fixed, we sum 
over all endpoints b such that bi < 62 < M. Hence one 
has 



P2,s{M,tm) = lim lim 



(51) 



^2,s(e,7?) 

where the probability weight is given by 
H^2,s(M-77,TM,e) = / d6 W^2(M-?7,TAf|6,e) 

J ord 

= / db / dy [5{y2 - {M - 77)) (52) 

J ord J Ol d 

X p<MAt>, TAf) P<Ma{y-, TAf |e, 0)| , 

where the integrations satisfy —00 < bi < h2 < M and 
~oo < yi < 1/2 < M ~ respectively. The normalization 
is such that 



/•l f'OO 

Z2,s{v,^) = / dTM / dAfM^2 



s(M-77,rM,e) ;. (53) 



As in the previous subsection, we will use only the leading 
order term of the expansion of W2.S and Z2^s in powers of 
r] and e. The computation of the probability weight can 
be done along the same lines as before (33 37), keeping in 



mind that only the final points change. At lowest order 
in r/, one has 



M 



62 00 

d6i -77^ 





dfc / dfc; 



-tm - 



(54) 



W2,s{M ~r],TM,e) = j db2 

— 00 

e-'^k'2k2^k,,k',{b)<i>kukMe-^''^''^^e 

where we have integrated over yi first, and then over k'l, 
with the delta function coming from the closure relation 
of eigenfunctions. For the eigenfunctions evaluated in the 
final points one obtains 



ki ,k: 



Xb) 



\/2 

(sin((7iCi)sin(g2C2) - sin(g2Ci) sin(giC2)) 

TT 



in terms of the previous qj = kjM, and with the variables 
1,2). Using the expansion in e 



= T?^ 



written in Eq. (38), one finds the leading order to be 
0{rf'() so that the expansion is 

W2,s{M-V:rM,e) = rj^e W2,s{M,tm) + 0{r,h,r,''e^) , 

(55) 

with 



W2,s(A//,rAf) 



r3M4 



dC2 



dCi 



dq 



<iq2 <i2q2^2{qi,q2)e 



X (sin((7iCi)sin(g2C2)-sin(g2Ci)sin(giC2)) , (56) 



where 62(91, 52) is given in Eq. (39). As before, the 



integration over the moments 52 and q2, corresponding 
to the top walker, can be factorized in the determinants. 
From the part [0,TAf], one recognizes T2(gi|Af, ta/), and 
for the part [tm, 1], we introduce 



f2{qi,C\M,l-TM) 



dq^ q^e 



X (sin((7iCi) sin(g^C2) - sin(g^Ci) sin(giC2)) • (57) 

With the help of the two functions T2 and T2, one ob- 
tains 



W2,s{M,tm) 
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7r3M4 



dC2 / dCi 



'0 JC2 

dqi T2{qi\M,TM)f2{qi,C\M,l - tm). (58) 



T2 and T2 can be written as determinants, as in Eq. (44 ), 
and 



T2(gi,C|M,t) = V^ 



M 
2t 



sin(giCi) 



sin((7iC2) 



e 2t 



(59) 



where we used the relation in Eq. ( E3 ) , and where t 
stands for 1 — ta/ (for compactness). The fact that we 
do not have the periodic boundary conditions introduces 
an asymmetry, thus the two determinants T2 and T2 do 
not have the same form. 

As before we also expand the normalization ^2,5(77, e) 
in powers of i] and e, the leading term being 

Z2,s{v, e) = '7^e ^2,3 + oiffe) , 

with K2^s a number, independent of M and tm which 
can be computed by the normalization condition: 



1 = 



/ dM / dTM P2,s{M,tm) 
Jo Jo 

= dM I dTMW2,s{M.TM) ■ (60) 

One finds K2^s — 2/-\/7ri and the joint probability distri- 
bution function, for two stars configuration is 



P2.,s[M,tm) 



it[tm{1-tm)?'^ 



tm „ ( M 

J^2 



UErff 

TM J V 



M 
71 



1 1 tm{1 - Tm) 

-e 2 



s/^M 2 - TM 
V2(^— ^1 ErflW^ — ^ I e'^<2-Af) 



l~TM M 



^_2-TMl fM 
A'P 2 \./2¥m 



(61) 
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FIG. 7: Distribution of the time at which the maximum is 
reached for N = 2 non-intersecting free Brownian bridges: 
the solid line is our analytic result given in Eq. ( 62 1 while 



the symbols are the results of our numerical simulation. The 
dotted line describes the pdf of the time tm for one Brownian 
path (whose cumulative distribution is given by Levy's arcsine 
law). 



Integrating over M one obtains the marginal, i.e. the 
pdf of the time to reach the maximum 



P2.s{tm)^ / dMP; 
Jo 



2,s{M,tm) 



arctan 



'tm(1 - Tm) 



2t 



M 



1 + TA/ V 1 - TM 



(62) 



which is plotted with a comparison with numerical sim- 
ulation in Fig. [7] One easily obtains the asymptotic be- 
havior of the distribution of the time tm'- 



• for Tm 0, one has 



P2,s{tm) 



3^2 1 



/2 



•m 



• for rj\/ — > 1 the probability distribution diverges as 



2 \/2 

P2,s{tm) ^ (1 - Tm) 



-1/2 



which can be compared to the divergence of the dis- 
tribution of the time tm for one Brownian motion 
(the Levy arcsine law, plotted in Fig.jT]) Pi^s{tai) = 

- [''■a/(1 — Tm)] ■ The fact that we have another 
uncrossing Brownian path below preserves the ex- 
ponent of the divergence, but changes slightly the 
pref actor. 

Notice finally the mean value is given by {tm)2.s — 7/2 — 
2y/2 = 0.671573 ... and that the sign of the derivative of 
P2,s{tm) changes in tm — 0.451175. 



FIG. 8: One realization of the N=2 excursions configuration, 
with the maximum M and the time to reach it tm- 



G. Periodic boundary condition and positivity 
constraint: N — 2 excursions configuration 

Now we come to the computation of the joint probabil- 
ity distribution of the maximum M and its position tm 
for the configuration of 2 vicious Brownian paths sub- 
jected to stay in the half line x > 0, in addition to the 
periodic condition: the two paths start in x — and end 
in a; = at r = 1. This corresponds to two Brownian 
excursions, with the non-crossing condition (we will use 
the subscript 'i?'). 

The positivity constraint is easily implemented in the 
path integral formulation by adding another hard wall 
in X = 0. The potential associated to each of the two 
fermions then reads 



cix) = 



if < a; < M 
+0O elsewhere . 



(63a) 



The eigenfunctions of the one-particle Hamiltonian 
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^Box(2^j)--^^+^Box(x,) 



are given by 

^nix) = 



^^^■sm^-x), n = l,2, 



with the associated eigenvalues 

1 n^n^ 

En 



2 M2 



(63b) 



(63c) 



(63d) 



This allows to express the corresponding propagator, for 
any point a, b and time < ta < ifc < 1, as: 



, (63e) 



Hl,n2>0 
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where n = (71,1,712) is a couple of two positive integers, 
and in which we use the Slater determinant ^'„(a) = 
deti<jj<2[V'«i(aj)]/y2!. 

Here the regularization is necesary for the boundary 
values of both fermions, because the lowest path cannot 
start and end exactly in a; = 0. Therefore we put a;(0) = 
a = (e, 2e) = e and a;(l) — b = e, where the notation 
£ — (e, 2e) differs slightly from the previous one e — (0, e). 
For this geometry, one uses Eqs. ( 30a|30c l to obtain the 
probability weight 

W2,e{M - ri, TM, e) = W2{M~ 77, tm\£, e) 
dy {S{y2 - {M - 77)) 

X PBox(e, l|y,TAf)pBox(y,TM|e, 0)| , (64a) 

where the ordered integral covers the domain < j/i < 
7/2 < M. The normalization function is 



^2,£;(»7,e) = Z2{'n\£,£) 



so that 



P2,e{M,tm) lim lim 



W2MM ~ ri,TM,e) 
Z2,e{v,<^) 



(64b) 



(64c) 



Inserting the spectral decomposition of the propaga- 
tor (63e) in Eq. (64a I, one has 



W^2,B(M-?7,TM,e) - / dy 5(2/2 - (M - 77)) 

J ord 



n.n' >0 



ord 

2,-. ^ , /2 



c*„(y)^';(£)e 2, 



■ (65) 



;s as 



We follow the same procedure as before in Eqs. p3||37 1 
using the symmetry in the exchange of indices (n-] o n: 
and n[ ^ n'2) we express the two Slater d 
the product of their diagonal components 

W2,£;(Af - r],TM,e) = 



[(l-TM)T.'^+XMrt^] 
2m2 



M-r 



X7/-„,(M-r,)7A„^(M-77) Jdyi t/.;, (7/i)7A„,(yi)}. 



(66) 

An expansion at lowest order in rj of this expression, mak- 
ing use of 



MM - 77) = (-1)"+!^/-— ,7 + 0{v^), (67) 



together with the closure relation of the eigenfunctions 

I'M 

/ dyi 7/;„,(7/i )?/';, (?7i) = (5„i,„'^, (68) 
^0 

yields 



W2,e{M -r^,TM,e) 



^ {*„.(£)*;(£) 



{-ir^+< 712 n'^ 6,,,,,^} . (69) 



The expansion to the lowest order in e of the eigenfunc- 
tions reads 

*r.(e) = (-1)^^ ^i"2 A2{nlnl) e* + Oie') , (70) 
where A2(Ai, A2) is the 2x2 Van der Monde determinant 



A2(Ai,A2 



1 Ai 
1 A2 



(A2-A1). 



Hence to lowest order, one has 

W^2,B(Af - 77, TM, e) =?7'e* W2,e{M,tm) + oirj^"") 
647rio 



W2,e{M,Tm) 



9M13 



EH 



(71) 



rii>0 



r!2>0 

X ^ (-l)"^7l2A2(772,772)e-5l^(l-"*^'"'| . 
ra3>0 

To ensure the joint pdf to exist, the normalization must 
have the same dominant term 

^2,B(?7,e) -?7'e' K2.E + 0(7^^^) 



Combining these expansions in Eq. ( 71 1 and in Eq. ( 72 ) 



one obtains from Eq. (64c 



P2MM,m)=^^ y: (-1) 



rii ,n2,n3>0 

^2 ^2^/„2 „2^„-;;^i^-[Tl?+rJ,./n^+(l-TM)ni] 



X (?^2 - Tlf )(773 - 77f )e 



(72) 



where the computation of normalization has been left in 
Appendix [Cj 

This formula is well suited for an integration over M 
to deduce the marginal law of tj\/: 



P2,e{tm)^ dM P2,e{M,tm) 
Jo 



1280 



E(-i 



K"2n3)^ {nl - nl){nl - nj) 

216 ■ 



ni>0 



K + '^Mnl + (1 - TM)nl] 



(73) 



14 



By difFerenciating N times the Jacobi identity 




FIG. 9: Distribution of the time at which the maximum is 
reached for N = 2 non-intersecting excursions (watermelon 
with a wall): the solid line corresponds to our analytic result 
while the symbols are the results of our numerical simulations. 



which is plotted in Fig. |9] 

Following the methods introduced in 19J, one is able 
to compute the asymptotic behavior of P2,e{t]\[) when 
tm — > 0. Starting from the expression of the joint dis- 



tribution of {M,tm) in Eq. (72 1, one writes the integral 
over M to get the marginal distribution P2.e{tm) as 

noo 

P2,e{tm)^ dM P2.e{M,tm), 
Jo 



and perform a change of variable 



2M2' 



(74) 



One obtains 



2^ 1 
StttI^ jo 



dx 

n2+n3^2„2^2/ 2 „2^^„2 „2^ 



P2.e{tm) 

ni, 712,713 

When Tm 0, the argument of the first exponential 
becomes —x{nl + n^)! \/tm at leading order, and the 
dominant term in the sums over n\ and ns is given by 
(ni — 1, = 2) and {n\ = 2, 713 = 1) because of the fac- 
tor (?i| — n^), which enforces n\ 7^ 71,3. Hence to leading 
order, one has 



2"^ 1 

P2.e{tm) 3 



dx x^ 



00 

E 

"2 = 1 



X V i-iy^nl {2nl ~ 5) e'^e 



(76) 



l+2^(-l)"e-"'- = 2(-) ^ ^e-'("+^)V^ (77) 

with respect to z, one has at leading order when z — )■ 

1/2 



(') (D "i.-- <™) 



With z = Xy/TM, one obtains at leading order 

P2ArM) - 2V/V;/^+^) /°°dx xV2,-(5^+fl)^^ 

Jo 

(79) 

which can be approximated by a saddle point method. 
At leading order when tm — ^ 0, one obtains 



^2,s(tm) 



Stt^ 1 

g-e • 

5 T^f 



(80) 



By symmetry one obtains the same behavior when tm — >■ 
1 by replacing tm by 1 — tm in the above formula. 

As in the previous paragraphs, it is possible to have a 
compact formula for the joint pdf, anticipating the gen- 
eral computation for N vicious Brownian paths. Let us 



introduce VI2, the analogue of T2 (431 in this discrete 
case: 

n2{n,\M,t) = Y,{-lfe^2{nl,e)e-^'''\ (81) 

with k a generic index standing for n2 (t ~ tm) or 713 
(t = 1 — Tm). Factorizing the sum over k in the last 
line of the Vandermonde determinant, and using sums 
defined by 



(M,t)-E(-1)'=A:2V 



ytk^ 



(82) 



fe>0 



one is able to write 

f^2K|M,t) 



uji{M,t) uj2iM,t) 



(83) 



Hence the joint pdf can be expressed in a factorized way, 
reminiscent of the formula for the watermelons case ( 42 ) 



X n2ini\M,TM)n2ini\M,l-TM). (84) 

In the following we will extend these formula, for bridges 
and excursions, to N particles, for generic N. 
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III. DISTRIBUTION OF THE POSITION OF 
THE MAXIMUM FOR N VICIOUS BROWNIAN 
PATHS 

In this section, we extend the previous method to 
any number N of vicious Brownian paths, ordered as 
Xi{t) < X2{t) < ■ ■ ■ < xn{t) for < r < 1 (again we 
take a unit time interval, without loss of generality). We 
want to compute the joint pdf of {M,tm), respectively 
the maximum and the time to reach it for different con- 
figurations, with 



M — max xn(t) 

0<T<1 



xn{tm)- 



(85) 



In the same way as in the previous section, the non- 
colliding condition translates into treating indistinguish- 
able fermions. They do not interact between each other 
but are all subjected to the same potential (which is one 
or two hard wall(s) potential in our computations). The 
Hamiltonian of the system is 



N 



(86) 



with identical one particle Hamiltonians for 1 < 

i < N. If ^Ei{xi) is the eigenfunction of ff^') associated 
to the energy Ei, one can write the eigenfunction of the 
iV-particles problem as a Slater determinant 

^e{x) = , ^det_^(VB,(a;j)) , (87) 



/J\[\ l<i,j<N 

with E ~ Ei and where we used the bold type for 
N vectors: x — {xi, X2, • . • , xn)- Hence the propagator 
reads as in formula ( 19 1, but with determinants of N x N 
matrices. The method to compute the joint probability 
density function is exaclty the same as before, replacing 
2-vectors by A^- vectors in formulas (30). 



A. TV vicious Brownian bridges 

We first consider the watermelon configuration, i.e. 
the case of N on-intersecting Brownian bridges (see defi- 
nition in subsection II E I . The periodic boundary condi- 



tions impose the starting and ending points to be at the 
origin. As discussed before, this needs a regularization 
scheme, and we take a = b = (0, e, 2e, . . . , {N — l)e) = e. 
The joint pdf is 

Wn^b{M - i],TM,e) 



lim 



where Wn,b{M — 1], tm, e) is the probability weight of all 
paths starting and ending in e with xn{t]\i) — M — rj, 
given by 



WN,B{M~V,rM,v)^ / dy\s{yN - {M 

Jord ^ 

X P<M,Af(e, l|y,TM) p<M,Jv(y,TA/|e,0) j. (89) 



P<M.N is the propagator of N indistinguishable fermions, 
each of them having the one-particle Hamiltonian 
H^j^i (29b), with eigenfunctions (j)k{x) and eigenvalues 



Ek given by Eqs (29d 29e) 



The starting point of our computation is given by 
Eq. (33) with the upper point is now y^v- AH the ma- 



nipulations done before for 2 paths extend to iV paths. 
A subtelty is the fact that the ordered integral over the 
domain —oo < yi < ?/2 < • • ■ < Un-i < Vn < M is pro- 
portional to the same unordered integral over the domain 
— oo<yi<Mforl<i<A^, because of the symmetry 
of the integrand. This is shown in Appendix [X] With the 
same argument of symmetry between the exchange of any 
ki -n- kj , and the same for k' one could use only the diag- 
onal terms for the two Slater determinants expressed in 
y. After permuting the y-integrals with the /c-integrals, 
this operation gives — 1 closure relations, i.e. a prod- 
uct of — 1 Dirac delta functions Hila^ H^i ~ '^Di ^iid 
the product 0fe„(M — rj)(f)i^i^{M — rf). The expansion of 
the Slater determinants evaluated at the initial and final 
points involves now the function 

^^('^^ " iif<N ('^^ ' ('^'^ ~ 4)) ' ^^^^ 

with the rescaled variables qi = kiM. After the afore- 
mentioned manipulations, one finds at lowest order 

WnMM - V, TM, e) = r;2e^(^-i) WN,BiM, tm) 

+ o , (91) 

together with 



Therefore one obtains from Eqs (91) and (92) together 
with Eq. (ISSl): 



Pw,b(M,tm) 



A 



N,B 



d(7i . . . dg7v_i e ~2m^ 



M^'+3 JO 

X TNi{q^}\M,TM)TNi{q^}\M, 1 - tm) , (93) 



where 

TN{{q^}\M,t) 



dqN qN e *^eAr(q) , (94) 



depends on iV — 1 variables {qi} = (gi, . . . ,qN-i)- In 
Eq. (93) the amplitude An,b, whose computation is left 



in Appendix [C] is given by 

22N-N/2 



A 



N,B 



(95) 



Moreover the function T n defined by Eq. ( 94 ) can be 
written as a determinant: 0jv is a determinant, and the 
integration over q^ can be absorbed in its last line. The 
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elements of this last line are indexed by the column in- 
dex j: 



dqN Qn^ cos (^qn - e 



/^AP+^\Jj{M,t) , (96) 



with the help of the family of functions 

Vj{t) EE \Jj{M, t) = t-^Hj ' 



(97) 



where Hj is the Hermite polynomial of degree j (see for- 
mula E2|. Then Tj^i{{qi}\M,t) is the determinant of a 



matrix, whose elements indexed by (a, b) are given by 



\ q^-^ cos (ga - &f ) if 1 < a < iV - 1 and 1 < 6 < iV, 
j y/^M''+'^\Jb{M, t) if a ^ and 1 < 6 < N. 

(98) 

The next step is to expand the determinant T^r in minors 
with respect to its last line: 

T7v(gi, ■ ■ ■,qN^i\M,t) = 

N 

^(-l)^+^'x/^AP+iU,(M,0 det [M^,,(ejv)] , (99) 
i=i 

where det [Mjv.j(OAr)] denotes the determinant of the mi- 
nor (Njj) of the matrix Qn, obtained by removing its 
N^^ line and its j*"^ column. These minors only involves 
the . . . , qN-i and M, not r^/- Expanding the product 
of the two functions T jv entering into the expression ( 93 ) 
yields 

TNi{q^}\M,TM) X TNi{q^}\M, 1 - Tm) = 

N 

(-l)'+^'7ril/r+^+2u,(TM)U,(l - tm) 

X det [MAr,i(eAr)] X dct [MaTj (Gat)] . (100) 

The dependence on the variables qi 's appears only in the 
last line of that expression (100), and the Cauchy-Binet 



formula (B7) allows to compute the integration with re- 
spect to 9i, . . . jQat-i of the product of determinants of 
minors: 



dqi . ..dqN^i e 



^k = l "k 



(101) 







X det [MAr^,(eAr)] det [MA,j(eAr)] 

= (iV- l)!det 



where a, b are respectively the line and column indices of 
the matrix of which we take the determinant, and Dn.b 
is the N X N matrix with elements 

(102) 



This matrix enters into the expression of the cumulative 
distribution of the maximal height of N watermelons [32] : 



Proba 



max xm(t) < M 

0<T<1 



1 



Putting together (100) and (101) in formula (93), and 



dividing by Kj^ b given in Eq. (92), one finds the joint 
pdf of M and tm 



det 'Djy^B- 
(103) 



N 



PnAM.tm) - Bn^b J2 (-l)'+'U,(rM)U,(l - tm) 
X det [M,,j (Djv,s)] , (104a) 



or equivalently 



Pn,b{M,tm) - -Bn,b det ( tu(i'^^'^^^) 



U(rM) 




Bn.b det [Bn.b] *U(1 - rM)D^isU(rAf) , (104b) 



with the normalization constant B^^^ — V^Yl^^i j!2-' . 

It is also interesting to characterize the small 
Tm behavior of the distribution of Pn,b{tm) — 
Jq PN,B{M,TM)dM . To study it, it is useful to start 
from Eq. ( 93 ) which, after integration over M yields 



^ 1 /"OO /"OO />oo 

Pn,b{tm) ^ An \\ / dqi / dqNqN / dq^q^ 
,_i Jo JO Jo 



6jv(gi,--- ,gjv)9(gi,--- ,g^) 

(1 - TMjq^ + TMqN^ + Lj=1 <lt 



, (105) 



with An = 2^'/2r (l + ^) An,b- To study the behav- 
ior of Pn,b{tm) when tm — >■ 0, we perform the changes 
of variables qi — ^TMVi, for 1 < z < iV so that when 
Tm ^ 0, one has to leading order 

' M III ^uiyN I dq'NqN 



Pn.b{tm)^t^j " 1 [ / dy.VN I dq'NqN 
,_i Jo Jo 



^Niy/TMyi, - ■ ■ ,^/rMyN)&{^/TMyi, - ■ ■ ,g^) 



I 2 , v^Af 2 



1+^ 



(106) 



Besides, one has, also to leading order 



0(r 



and similarly 



M 



(107) 



+ 0{t^j - +') . (108) 
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with some regular function f{q'pj). However, being anti- 
symmetric under the permutation of the variables j/^'s it 
is easy to see that this leading order term ( 107 108 jvan- 
ishes after the integration over yi, - • ■ ,yN in Eq. (106 1. 



In fact if one writes the low tm expansion of the product 

as 



^^r(^/T5v^2;l, 

jV^ + (jV-l)^ 
— T ^ 



[Vq + VlTM 



(109) 



where z/j,'s are functions of yi, ■ ■ -yNjlN '^^^ '^^^ show 
that the first non-vanishing term, after integration over 
these variables, corresponds to Vk with fc = 1 for = 2, 
fc = 2 for = 3 and more generally fc = iV — 1 for generic 
iV. Therefore one has, to leading order 



Pn,b{tm) 



' M 



(110) 



although the explicit calculation of the amplitude is a 
very hard task. 



B. A'^ vicious Brownian excursions 

We now consider watermelons with a wall, i.e. N ex- 
cursions under the non-crossing condition (see the defi- 
nition in paragraph |II G I . Again the starting and ending 
points arc at the origin, but with this configuration we 
must use a = b = (e, 2e, . . . , Ne) — e, which is slightly 
different from the cases of bridges, because now, all paths 
must be positive. The joint pdf of M and tm is 



Pn.e{M,tm) = lim 



(111) 



where the probability weight of all paths starting and 
ending in e, verifying for all < t < 1 the non-crossing 
condition < a:i(T) < • • • < xn{t) < M, and for which 

xn{tm) ^ yN ^ M - r], reads 

WN,EiM -T],TM,e) ^ / dy\5{yN -{M ~rj)) 

Joid ^ 

X PBox(£, l|y,TM) PBox(y,TAf|£,0)|. (112) 

Here the integration is over the ordered domain < j/i < 
• • • < yN < M. The boundary conditions of the ex- 
cursions configuration enforce the use of the propagator 
PBox associated to N fermions confined in the segment 
[0,M] (63e): the Hamiltonian is of the form (86 1, with 



the one particle Hamiltonian ( 63b ) . 

One can follow the first steps of the paragraph |II G 



replacing now 2-vectors by iV-vectors. As in the bridge 
case, the ordered integration can be set unordered {i.e. 
integrating over < < Af — 77 for all 1 < i < — 1), 
with a numerical prefactor (see appendix [A|) . The spec- 
tral decomposition is done with sums over 2A^ indices 



, n^v) and n' = {n[, 



(■TV 



), as in for- 



mula (65 1. Using the symmetry of the integrand under 



the exchanges of any couple Hi -O- rij and n- O n'j (for 
1 ^ ijj < one can write the Slater determinants of 
the intermediate point (at t = tm) under their diago- 
nal form, with a numerical prefactor, irrelevant in the 
computation. Hence one has the formula 



Wn,e{M - i],TM,e) 

JV(JV-l) 



A^!2" 



E 



e 2 M ■ 



r [(l-rM)Ti'^+rA/7 



X ^-^^ (£)*„(£) 7/^„^(M-77)^„„(M-77) 
N-l C 

/ dyj Vn;(yi)'0n,(2/») )• (113) 



n 



where ipn{x) is given in Eq. (63c). The expansion to 
lowest order in 77 gives A^ — 1 closure relations (the last 
line of the preceding expression) 



<^y^ '4'n'Xyi)i^ni{yi) = 5„ 



(114) 



and we use formula ( 67 1 to get 



Wn,e{M - ■r],TM,e) 



^ (-l)"~+"'«„^„^ 



N-l 

■ E ' 

i=l 



e~5i72 [(1— ^m)"W^+-^m"n. 
X det?/'Ti(e) det^fi(e) 



(115) 



where h = (ni,n2,-- - ,n'j^). One can then use the ex- 
pansion to lowest order in e 



detVr.(e) = (-1)" 



N 

n 



2^7r 



l[n, A^^(n?,...,n^) +0(6^^), 

(116) 

where AAr(Ai, . . . , Aw) is the Nx N Vandermonde deter- 
minant. In this formula the first line contains only nu- 
merical constants, irrelevant for the computation (they 
will be absorbed in the normalization constant), only the 
second line is relevant, containing the dependence in M 
and in the indices of the sums Ui, and the lowest order 
in . Hence one has to lowest order 

WnAM - V. TM,e) = £2^' Wjv,b(M, tm) 

+ 0(77V^',77V^'+i) , (117) 

as well as, for the normalization, 

Zn.Mv, e) = v'e'"' K^,e + o(^'e'^') , (118) 
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with Kn^e a number. Taking the double hmit e, — ^ 
one obtains the joint probabihty 



PN,EiM,TM) 



'^n.e{M,tm) 



K 



N.E 



A 



N.E 



N-1 
2 TT 2 



]ljN{2N+l)+3 



X Ajv(ni, . . . ,n^_i,n^) AN{nl, . . . ,n%_i,n'j^ ) 

e~^[(^"'''-f'"'"'+^^f""] I , (119) 



iV-l 



X e '=1 



where the numerical constant is given by (see ap- 
pendix C 2 1 



A 



N,E 



(120) 



Following the notations introduced in the N ~ 2 case, 
we define [writing for compactness the A^— 1-uplet {n-i} = 
(ni, . . .^npf-i)] 



npan,}\M, t) = £(-i)'=fc2A^({„2}, 



rtk^ 



k=l 



1 nf 
1 nl 

U>l{t) UJ2it) 



2(JV-1) 
'l 

2(JV-1) 



WAr(t) 



(121) 



where we used the property of the Vandermonde deter- 
minant to factorize the sum with respect to k in the last 
line of the determinant, which plays the role of as- 
sociated to t = tm, or the role of with < = 1 — tm, 
resulting in elements like oJi {t) = u!i {M, t) , for 1 < i < 



(122) 



fc=i 



With this notation, we obtain the compact formula 



Pn^e{M,tm) 



A 



N.E 



]]^fN(2N+l)+3 



HE 



i—l \_ni — l ) 
X r!w(K}|M,TM)f^Ar({n,;}|Af, 1 - tm) , (123) 



which is the "discrete" analogue of formula ( 93 ) obtained 
in the case of N bridges. 

The ultimate step is to expand each determinant rj^r 
with respect to their last line. Then, using the discrete 
Cauchy-Binet identity with respect to the sums over n^, 
one has 



Pnm{M,tm) 



B 



N 



M2 \KP 



N,E 



X a;,(TM)wj (l - tm) det [Mij(Djv,B)] , (124) 



with the normalization 

(-1)^+177^/2^2 2-2^' 



B 



N,E 



l+j) 



(125) 



and where appear the minors of the matrix Dn,e 
^n,e{M) whose elements are, for 1 < i, j < N 



+ 00 



^^^E^,,= E i?2(.+,-i)(V2Afn)e-2*^'" . (126) 



This matrix enters into the expression of the cumulative 
probability of the maximum |29| : 



Proba 



max xn(t) < M 

0<T<1 



= (-1) 



N 



det Dn E 



nL(2j-i)! 



Defining the vector with elements, for 1 < i < N 

27r2 



VE^{t) = \JE^iM,t) 



1 

M 



AP 



(127) 

Lj,{t) , (128) 



one is able to express the result as the determinant 



PN,EiM,TM) 



Bm,e det 



^N,E ^e{tm) 



*Ub(1 -tm) 







or as the matrix product 



(129a) 



Pn,e{M,tm) 

= Bn,e det [D^v,b] 'Ve{1 - TM)r>Jj'E^ e{tm) ■ (129b) 

The behavior of the marginal distribution Pn.e{tm) 
when Tm — >■ can be obtained as explained in para- 
graph |II G| for = 2 vicious excursions. The leading 
behavior, when tm 0, is given by 



Pn,e{^m) 



^4 

Am V 



JV(iV + l)(2W + l) 



(130) 



However the algebraic correction to this leading behavior 
is hard to evaluate. 



IV. NUMERICAL SIMULATIONS 

In this section we present the results of our numerical 
simulations of the distribution Pn,b/e(j'm) both for the 
cases of bridges and excursions. To generate numerically 
such watermelons configurations we exploit the connec- 
tion between these vicious walkers problems and Dyson's 
Brownian motion, which we first recall. 



A. Relationship with Dyson's Brownian motion 

To make the connection between the vicious walkers 
problem and Dyson's Brownian motion, we consider the 
propagator of the N vicious walkers 

'Pjv(v,t2|u, ii) = VNivi, ■ ■ ■ ^VN^hlui, ■ ■ ■ ,UN,ti) , 

(131) 
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which is the probability that the process reaches the con- 
figuration 2:1(^2) = wi, • • • ,xpf{t2) = vn at time t2 given 
that Xi{ti) — Ui,--- ,xjv(ii) — uj^ at time ti: this is 
thus a conditional probability. This propagator satisfies 
the Fokker-Planck equation: 



'Pjv(v,i2|u,ti 



1 ^ 



1=1 



d 

dv. 



2 7'Ar(v,t2|u,ti) , (132) 



together with the initial condition 

^Ar(v,t2 = <l|u,<i) = 5W(v-u) 



(133) 



and the non-crossing condition, which is specific to this 
problem, 

7'jv(v,t2|u,ii) = 0, siv, =vj ,V(ti,t2) • (134) 

On the other hand, let us focus on Dyson's Brownian 
motion. For this purpose we consider random matrices 
whose elements are time dependent and are themselves 
Brownian motions. For instance, for random matrices 
from GUE, i.e. with /? = 2, we consider random Hermi- 
tian matrices H = H{t), of size N x N, which elements 
are given by 



^ (brnnit) + 
Hrnn (0 \ ^mm (^) i 



71 ^ 



.it) 



.it) 



m < n , 

m — n (135) 

m > n , 



where b„in{t) Si,nd 6m„(i) are independent Brownian mo- 
tions (with a diffusion constant D = 1/2). We denote 
Ai(t) < A2 (<)<■■•< AAr(t) the N eigenvalues of H{t) 
{ 135[ ) [or more generally of a matrix belonging to a /3 
ensemble with /3 = 1,2,4 which is constructed as in 
Eq. ( 135[ ) with the appropriate symmetry]. One can then 
show that the A^'s obey the following equations of motion 
(which define Dyson's Brownian motion) [41] |42] 



dX^jt) 
dt 



^ Kit) 



Kit) 



mit) , (136) 



where rjiS are independent Gaussian white noises, 

{V^it)^l3it')) - 5^j5it ~ t'). Let 7'Dyson(A, ^1/^, 0) = 

7'Dyson(Ai, • • • , Aat, • • • , ^XN, t = 0) be the propaga- 
teur of this Dyson's Brownian motion ( |136 ). It satisfies 
the Fokker-Planck equation 



d_ 
dt 



N 



2 ^ aAf 

Z— 1 ^ 



-It' 



A 



yson 



(137) 



Of course, one has also here 7^Dyson('^, t\^, 0) = if A; — 
Xj , for any time t and from Eq. ( 137 1 one actually obtains 



^Dyson(A,t|/X,0) - (Aj - A,)^ , A, ^ Xj 



(138) 



This Fokker-Planck equation can be transformed into a 
Schrodinger equation by applying the standard transfor- 
mation [751: 



7'Dyson(A,<|/X, 0) = 



exp 


2 ^l<i<j<N 


(A, 


-A,) 


exp 


2 J2l<i<j<N 




- 



X WDyso„(A,t|M,0) , (139) 
where VVDyson(A, t\fi, 0) is such that 

WDy«o„(A, t = 0|/x, 0) = ,5(^) (A - m) . (140) 



On the other hand, given Eqs ( 138 1 and ( 139 1 one has 



WDyso„(A, t\fl, 0) - (A, - A,)'3/2 ^ ^ (^4^) 



From Eqs (137) and ( |139 1, one obtains that WDyson sat- 
isfies the following Schrodinger equation 



yson — r, 2^ 



Ij^WDyson - ^ 



2 ^ aA 



2 WDyson 



TV 



(A, - AO^ 



WDy 



.(142) 



Hence for a generic value of /3, ( 142 ) is the Schrodinger 
equation associated to a Calogero-Sutherland Hamilto- 
nian (on an infinite line) [77l |78] . However for the spe- 
cial value (3 = 2, the strength of the interaction van- 
ishes exactly and in that case the equation satisfied by 
WDyson ( 138 140 142 1 is identical to the one satisfied 



by the propagator in the vicious walkers problem ( 132 
133 134 1 : in that case Eq. (142) corresponds to the 
Schrodinger equation for free fermions. In other words, 
for the special value /3 = 2 there exists a simple rela- 
tion between vicious walkers and Dyson Brownian mo- 
tion which can be simply written as 



'PDyson(A,t|/X, 0) 



n.<,-(Aj--Ao 



VNi\t\f^,0). (143) 



In particular, for watermelons without wall on the unit 
time interval, for which the initial and final positions 
coincide, this relation (143) tells us that this process 



is identical to Dyson's Brownian motion where the el- 
ements of the matrices are themselves B rown ian bridges 



i.e 
B 



— > Bij in Eq. (1351 such that 



„(0) = B,,il) = and B,j(0) = 5,^(1) = 0. 



B.jil) = 0. This then 
allows to generate easily watermelons because one can 
generate a Brownian bridge Bijit) on the interval [0,1] 
from a standard Brownian motion bij it) via the relation 
Buit) = h,it)-th,il). 

For non- intersecting excursions, one can show that the 
dynamics corresponds to Dyson's Brownian motion asso- 
ciated to random symplectic and Hermitian matrices, as 
explained in Ref. [25] . 
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B. Numerical results 



To sample non-intersecting Brownian bridges we use 
the aforementioned equivalence to Dyson's Brownian mo- 
tion. We consider the discrete-time version of the matrix 
in Eq. (fT35|, at step < fc < T: 



Bu{k) 



1 

72 



I <i < j < N, 
i=3,l<i<N, 
I <i < j < N, 

(144) 

where Bij{k) are independent discrete-time Brownian 
bridges. They are constructed from ordinary discrete- 
time and continuous space random walks bij{k) where 
each jump is drawn from a Gaussian distribution. The 
aforementioned relation, to construct Brownian bridges, 
reads in discrete time Bij{k) — hij{k) — {k/T)hij{T). 

Diagonalizing the matrix Hij{k) at each step k gives 
the ordered ensemble of N distinct eigenvalues Ai(fc) < 
A2(fc) < ••• < Xjf^k), which is equivalent to a sam- 
ple of non-intersecting Brownian bridges, see Eq. |143| 
For each sample, we keep in memory the maximal value 
of Xpf(k), and the step time km at which XNikivi) = 
maxo<fc<T AAr(fc). The plots shown here are obtained by 
making the histogram of the rescaled time km /T (with a 
bin of size 8 to smoothen the curve), with T = 256 steps 
averaged over 10^ samples, for N = 2,3,... ,16. For 
Brownian excursions, one has to use another symmetry 
for the matrix, which is explained in Ref. |29j . 

First we present the comparison between our exact 
result for Brownian bridges and numerical simulations, 
for TV = 2,3,4,5,10. The agreement between our nu- 
merical results and our exact analytical formula is per- 
fect, see Fig. [TOj Then, anticipating on the applocation 
of our results to growing interfaces, we show a plot of 
the rescaled distribution as a function of the rescaled 
time um = {tm — l/2)/N^^/^ . Despite the fact that 
this scaling requires in principle that 3> 1, one sees 
that the rescaled distributions for moderate values of 
N = 10, 12, 14, 16 fall on a single master curve, see 
Fig. |11[ However, given the relatively narrow range of 
N explored in these simulations, one can not be sure 
that the asymptotic, large A^, behavior is reached. 

To compare our analytical results for A'^-watermelons 
to growing interfaces models and the DPRM, we have 
simulated both the PNG model and a discrete model of 
DPRM on a lattice. To sample the PNG model, we have 
applied the rules that were discussed before (see Fig. [T]) 
where the time interval is discretized. In particular we 
use a discretization of the nucleation process: at each 
time step, a nucleation occurs at site x such that \x\ < t 
with a small probability p, whose value is chosen such 
that one recovers a uniform density of nucleation events 
p = 2 in the continuum limit. The interface evolves dur- 
ing a time T . We compute the distribution of the position 
Xm of the maximal height of the dro plet (see Fig. [2]left), 
averaged over 10^ samples. In Fig. 12 we have plotted 
the distribution of the rescaled position um = Xm/T'^^^ 
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FIG. 10; Plots of the distribution Pn.b{tm) of the time tm 
to reach the maximum for A'^ = 2, 3, 4, 5, 10 non-intersecting 
Brownian bridges. The dots correspond to our numerical data 
while the solid lines correspond to our analytical predictions. 
There is no fitting parameter. 
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FIG. 11: Scaling form /jv, s (mm) of the distribution P]v,B(rM) 
with TV = 10, 12, 14, 16. The x-axis is the rescaled time um = 
{tm — |) X N^''^, and the y-axis is the rescaled distribution 



fN.B{uM) = iV"'/^ X Pn,b{1+umN- 



for T = 64 and T = 90 (both of them being statistically 
independent): the collapse of these data on a single mas- 
ter curve is consistent with the expected KPZ scaling in 

On the other hand, we have simulated a directed poly- 
mer model in 14-1 dimensions as depicted in Fig. [3] (left), 
where on each site, there is a random energy variable, 
with Gaussian fluctuations. We find the optimal path 
of length L , which is the polymer with minimal energy, 
with one free end-point. And we compute the distribu- 
tion of the position Xm of this end-point, averaged over 
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10^ samples. In Fig. 



12 we have plotted the distribution 
XAf/(aL2/3) for L 



of the rescaled position um = XM/{aL^' ) for L = 100 
and L = 200 (both of them being statistically indepen- 
dent): the collapse of these data, for different L, on a 
single master curve is consistent with the expected KPZ 
scaling in L'^/'^, while the value of a ~ V2 was adjusted 
such that these data coincide with the ones for the PNG. 

For comparison we have plotted on this same graph, 
on Fig. [12] the data for the rescaled time um of the 
maximum for N non-intersecting Brownian bridges, i.e. 
the same data as in Fig. 



11 



(except N — 10), without 
any additional fitting parameter. One sees that the fi- 
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Ui,: rescaled time ii, or position X^, 



tegral approach for free fermions, which incorporates in 
a rather physical way the non-colliding condition. We 
have considered three different types of vicious walk- 
ers configurations: non-intersecting free Brownian mo- 
tions, "stars" configuration, for N — 2, non-intersecting 
Brownian bridges, "watermelons" configurations, for any 
number N of vicious walkers and finally non-intersecting 
Brownian excursions, "watermelons with a wall" also for 
any N. We have also checked our results using numerical 
simulations of Dyson's Brownian, which uses the deep 
connection between vicious walkers and Random Matrix 
Theory, which we have reminded in detail. 

Thanks to the connection between vicious walkers and 
the Airy process ([t]), our results yields, in the large N 
limit, results for the distribution of Xm, the position of 
the maximal height of a curved growing interface in the 
KPZ universality class (Fig. [2] left) or the transverse co- 
ordinate of the end-point of the optimal directed polymer 
(Fig. [2] right). We have shown that our analytical results 
for finite N ^ 15, correctly rescaled, are in good agree- 
ment with our numerical data for the PNG model and 
DPRM on the square lattice. Performing the large N 
asymptotic analysis of our formula (93 123) remains a 



formidable challenge, which will hopefully motivate fur- 
ther study of this problem. 



Acknowledgments 



FIG. 12: Distribution of the rescaled time (respectively po- 
sition) Um for N vicious Bridges (respectively for PNG and 
DPRM models, see Fig. [2| on a log-lin plot. The dashed lines 
correspond to the PNG model, the symbols to the DPRM, the 
solid lines to the watermelons and the dotted line is a Gaus- 
sian fit. See the text for more details. 



nite N results describe quite accurately the center of the 
distribution, whereas there is a systematic deviation in 
the tails of the distribution (as expected since N is rel- 
atively small). Both data, for watermelons on the one 
hand and for the PNG on the other hand, suffer respec- 
tively from finite TV and finite T effects. Our data for 
the PNG seem to indicate that when T increases, the 
probability weight is slightly displaced from the tails to 
the center of the (rescaled) distribution: this suggests 
that the limiting distribution is not Gaussian (the best 
Gaussian 



12) 



fit is plotted with a dotted line in Fig 
This convergence to the asymptotic distribution is how- 
ever quite slow. 



V. CONCLUSION 

In conclusion, we have presented a detailed derivation 
of the analytic computation of the distribution of the 
time at which the maximal height of N vicious walkers 
is reached. Our analytic approach is based on a path in- 
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Appendix A: Integration over ordered variables 

In the main text, we have to sum over realizations with 
the constraint xi < X2 < ■ ■ ■ < xn~i < xn- Here we 
study how integrals may be simplified by the symmetry 
of the integrand. 



1. N — 2, two variables 

Consider a box [a, b] in one dimension, and a function 
of two variables f{xi,X2). In the main text, we want to 
compute ordered integrals like 



/ dX2 / dXi/(xi,X2) 
J a J a 



(Al) 



This is usually easier to compute the integral over the 
whole domain, that is to say 

/ dx2 / Axif{xi,X2) ■ (A2) 

J a J a 

Here we make a link between those two quantities, when 
/ is symmetric in the exchange of Xi and X2 f{x2,xi) = 
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f{xi,X2)- Let us start from equation (A2) with a sym- 
metric function / 



/ dX2 / dxif{xi,X2) 
J a J a 



b PX2 rb nb 

dx2 / dxif{xi,X2)+ / dx2 / da;i/(a;i,a;2) 

a J a J a J 

rb PX2 rb rxi 

/ dx2 / da;i/(a;i,a;2) + / dxi / da;2/(a;i, 0:2) 

'J a J a J a <J a 

PX2 rb nxi 

dx2 / dxif{xi,X2) / dxi / dx2f{x2,xi) 

<J a 'J a J a 

=2 / dX2 ( dxif{xi,X2) (A3) 
J a J a 

From the second to the third hne we use triangular inte- 
gration {a < X2 < b, X2 < xi < b} = {a < xi < b, a < 
X2 < Xi} in the second term. From the third hne to the 
forth, we relabehed the variables in the second term. In 
this operation if / was antisymmetric, we would have a 
minus sign, and hence a cancellation of the two terms. 
The conclusion for two variables is that if / is antisym- 
metric, the integration over the whole domain is zero, but 
if / is symmetric in the exchange of the two variables, we 
have the useful formula 

pb /•X2 ^ rb rb 

/ da;2 / dxif{xi,X2) = - / dx2 / da;i/(xi, a;2). 

'J a J a <J a J a 

(A4) 

2. A'^ variables 

The previous result can be extended for N variables, 
and the result is, for /(I'l, . . . , x^r) symmetric in the ex- 
change of any of two of its arguments 



/ dxfix) 2 — / dxfix) , 

J Ol d J 



(A5) 



with X — {xi, . . . , xn) and an integration on the domain 
a < xi < X2 < ■ ■ ■ < xn < b in the left-hand side (Ihs), 
and an integration over the whole domain a < Xi < b for 
1 < i < TV in the right-hand side (rhs). 

The proof is done by applying iteratively 
the formula ( |A4[ ) first to the — 1 couples 
{xn,Xn-i), {xn,xn-2), .■.,{xn,xi) changing the 
upper limit of integrations from b to xn ■ This operation 
results gives a factor 2 for each couple so that: 



b pb j*b pb 

dxN / dxjv-i... / da:2 / dxi f{x 

a J a J a J a 

b rXM fXN fXjv 

dxN da;Ar_i 



=2 



Af-1 



a 'J a 



rXN nXN 

/ dx2 / dxi f{x) . 

J a J a 

(A6) 



Then one repeats this procedure for the N — 2 cou- 
ples {xn^i,xn-2) to {xn-i,xi), which gives a factor 



2^ ^, and continues until step k with the N — k cou- 
ples {xN-k+i,XN-k) to {xN-k+i,xi) to have 



pb pb 




/ dxN / 


da;Ar_i 


J a J a 




2N—I2N — 2 


2N-k 




i-xn 


X dxN 


/ dx 


J a J 


'a 


rXN-k + l 






dxN-k- 


J a 





dx 



ordered 



dxi fix) (A7) 



unordered 



Doing so until fc = A^ — 1 one obtains the result ( A5 ) 



Appendix B: Cauchy-Binet identity 

1. Continuous version 

In the main text, we are confronted with integrals of 
two determinants that takes the form: 

jW=/da; det_ det_(gk{xi)), (Bl) 



l<i,i<JV 



l<k,l<N 



where da; = YliLi 1 5: * !i A^ are functions of 

one variable, and each variable Xi has the same domain 
of integration (in the main text we usually integrate over 
the momenta ki from to -l-oo). We have 



N 



N 



N 



i=l 



N 



o-'eSw fc=i 

N N 



ctGSjv i=l i=l 



N 



Y Y[9kiXa'{k)) 
.cr'GSjv 



fc=l 



(B2) 



Sn is the symmetric group, e(cr) is the signature of the 
permutation a e Sn- Now we make the change of vari- 
ables x^ = a;o-(i)- That is Xj = ^^^-iQ) ^^'^ Jacobian 
is unity. So we have 



N N 



/W= Y eia) fUdx^llMx'^ 

N 



k=l 



(B3) 
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is a group, so it exists a" € Sn so that a" = 
a^~^^ o a', and £{a") = e{a~^)e{a') = e{a)e{a'). More- following relation, valid for any integrable function h{x 



Moreover, one can easily deduce from Eq. (B7| the 



over Eaes^ = E.-io.es„ = E^"es«- Getting back to 
the Xj variables we have 



N N 



N 

t'^Sn i—k 
» N 

N 

. W N 



k=l 



E '^i^")Yl9k{x^„(^k)) 



L<T"eSA 



fc=i 



(B4) 



In the last line, the term in bracket is simply the deter- 
minant deti<k.i<Nigk{xi)). The permutation a does not 
appear anymore, so we have X^o-eSjv ^ ~ CIard(5'iv) — N\. 
Hence 



. N N 

Nl / nd2^*n/'(^^) ,<if<j,^9k{xi)). (B5) 
i=l i=l - ' - 



At this stage, the basic result is that we have expressed 
the first determinant only as a product of the diagonal 
terms. This is used in the main text. 

To go further we use the fact that the determinant of 
one matrice is the same as the determinant of its trans- 
pose: 



N N 



- TV! f fldx,f[h{xk) 

^ i — 1 / — 1 



1=1 k=\ 



N 



E ^'^'^)^9<y{k)[Xk) 



N 



N 



k=l 



= N\J2 J tl<ix^f[{fk{xk)gaik){xk)) 

(tGSn "i— 1 A: — 1 

= ^! E ^('^) n I / ^^fki^)9.ik)ix)\ ■ (B6) 

(TSSiv fc=l ^ ^ 

From the second to the third line, we factorized the in- 
tegrals, and we wrote x — Xk for each k. We recognize 
the determinant of the matrix whose elements indexed 
by are given by an integral of the product fi x gj^ 
yielding the so-called Cauchy-Binet identity: 



Ax det ifiixj)) det {gk{xi)) 

l<i,j<N ■' l<k,l<N 



^^■i<t<N^ I ^^M^)9jix)}. (B7) 



/ da; TT/i(a;,) det Mxj)) ( det g,{xj)) 



m ^^det^^ <{ / dx hix)Mx)g,{x) } . (B8) 



This relation is used in the main text. 



2. Discrete version 



Here we consider the following multiple sum: 

= E E • • • E i<^f<^ (M^j)) i<drt ^ (g,(n,)) 



rii n2 njv 



N 



^E<-)EE---En/'K«) 

^ ,^det „(g„.(a;j )) . 



(B9) 



Changing the indices of the sums mi — ^^(i), or rik — 



N 



in=j:<-) e e ■■• e n/»(-^ 



TV 



= E<-)EE---En/-(-^) 



l<zj,<A'^ 



(BIO) 



where we re-arranged the columns in the remaining de- 
terminant in the 'natural' order, thus factorizing the sig- 
nature of the permutation a. Hence there is a sum over 
the permutations a G Sn of £^(0"), giving the cardinal A^!. 
We obtain as in the continuous case that the first deter- 
minant is written as the product of its diagonal terms 
with a factor iV!: 

N 

= ^! E E • • • E n /'(-^) i<df<^ 

mi m2 mjv i—1 

N N 

=^!EE---En/™.(^^) E <^)n5^w("^^) 

mi m2 mjv «— 1 ctGSn "i— 1 

TV 

= N\ J2 <^)EE---En(/'(^"^)5^w(™')) • 

(Bll) 
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Factorizing the sums, one gets 



= iV! det (^/.(m)g,(m)) . (B12) 

\ m / 

In the same way as in the continuous case, if one considers 



where the amphtude aN,B can be computed using a Sel- 
berg integral [30] 



■)2N 



(2^)^/^nf=iJ 



(C3) 



The pdf of the maximum, F'j^ g{M) = dMFM,B{M) is 
thus given by 



^ .J'^iM .J'^iM (-9^^"^)) ' (^13) 

l<i,j,<N l<i,j,<N 

one finahy obtains 



Jn = N\ det ( V h{m)f,{m)gj{m) | . (B14) 

l<i,j<N \ ^ — ' / 



Appendix C: Normalization constant 

In this section, we compute the normaUzation constant 
of the joint distributions Pn,b{M, tm) and Pn,e{M, tm) 
for non-interesecting bridges and excursions respectively. 

1. Non-intersecting Brownian bridges 



To compute the amplitude s in Eq. ( 95 1 we use 



the fact that the integration of Pn,b{M, tm) m Eq. (93 1, 
i.e. Jq Pn,b{M, tm)<^tm must yield back the pdf of M 
which was computed in Ref. [SU]. Indeed in Ref. [30] it 
was shown that 



Fn,b{M) = Proba 



max xm(t) < M 

0<r<l 



(CI) 



dqe-^ [QNiq)f , (C2) 



FkBim 



-N 



2 0-N,B 



dq e 



X cos (gp - det[M,- ^(ejv)] det[M,,w(ew)] , (C5) 
where, in the last two lines, we have expanded both de- 
terminants in minors with respect to the last column (we 
recall that Mjjv(Ojv) denotes the minor {j,N) of the 
matrix Qn, obtained by removing its j*'' line and its 
N^^ column). Note that these minors do not depend on 
qn- This suggests to consider, in the last two lines of 
Eq. (C4), only the integral over the variable qj\j where 



we perform an integration by part. This yields: 



QN -In / .7r\ / 7r\ 

9^ cos (^^AT - j-J cos (^gw - i-J 

(j + ^ - 1) j d^AT cos [qN - j-j 



xcosl^qN -l-je + j dg^v^jv ^ 21/2 







xd, 



Qn 



(gjv - j|) COS (^qN - 



(C6) 



Inserting this result of the integration by part ( C6 ) in 



Eq. ( C4 ) and using the following identity (which we have 
explicitly checked for N — 2, 3, 4, 5 but can only conjec- 
ture for any integer > 5): 



Jo Jo 



sin (qi) 
qi cos (gi) 
sin (gi) 



gf-^cos (gi-iVf) 

sin (gi) 
gi cos (gi) 
9i sin (gi) 



sin (gAf) 
gjv cos (gAr) 
g^sin 



(in) 



q^-^cos (gjv-iVf) 

sin (gAr) 
2gAr cos (gAf) 
3g^ cos (gAf) 



gf-icos(gi-7Vf) ••• iV gj^-^ cos (g^ - iVf ) 
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one obtains finally 

-NaN,_ 



N 



This integral (C14) can be easily evaluated using 
residues, this yields 



j.i=i 



gj{s)^-{M^2sy-\ 



(C15) 



det[M/_Jv(^'Jv)]sin [2qjv - {I + j)-) ■ Using this result ( |C15[ ) in the expression above ( |C13[ ), 

^ 2 / nnr* nKf dine 

(C7) 

N-l ^oo 2 N 



On the other hand, from the definition of Pn^b{M, tm) 
given in Eq. ( 93 ) one has 



one obtains 
P(M,s) 



PN,B{M,TM)dTM = Fi,^siM) 



(C8) 



det[Mj-Ar(6'Ar)]det[M,,A,(6'N)](M\/2i)^+'-2e-2Mv^ 



(C16) 



with 

Pn,b{M,tm 



A 



N.B 



dqi... dqN-i e 



X TN{{q^}\M, TM)rN{{q^}\M, 1 ~ Tm) ■ (C9) 



Using the results identity shown above (C14 C15) we 
can write P(M, s) as 

i=i 



To compute this integral over tm ( C8 ) we compute in- 
stead 



JV 



A 



X Tjv ({<Z» } I M, TM)rN{{q^}\M,T- TM 

and clearly 



dgi . . . d(7jv_i e 



E (-lF+'det[M,-Ar(0^)]det[Mz,Ar(0Ar)] 



X 



(CIO) 



dgg- 



^_^.^,_, sin(2g-(j+/)f; 



(C17) 







2M2 



Under this form (C17) it is easy to invert the Laplace 



J PN,B{M,TM)dTM^P{M,T^l). (Cll) r 

. . Jo 

The structure of P(M, T) in Eq. (CIO) suggests to com- 



transform to obtain P(Af, T) and finally one obtains 
1 





PN.B{M,TM)dTM 



pute its Laplace transform with respect to T 



N 



P{M,s)= I P{M,T)e-'^dT 



3,1 = 1 

TT 



xdet[Mj,jv(6'jv)]det[Mi,Ar(6'jv)]sin 2gAr - [l + j 



A 



N,B 



N-l „oo 



lim , ^ 



n / 



-«l|9Jv|-e2kNl 



(C18) 



X / dgAT / dqNqNQN 



,0Af(gi, • ■ ■ ,9^) 



0Ar(9l, • • • ,qN) 



2M2 



(C12) 



Finally, usin g the iden tity ( C8 1 together with the explicit 

(C19) 



expressions (C7 C18) one obtains 

I 22N-N/2 



B — —Naj\i g — , 



where we have regularized the integrals over q^ and q'^ 
with the help of this exponential term e"'^"^''"'"'^^!''"! . 
The next step is to expand the determinants in Eq. ( C12 ) 
with respect to their last column. This yields 



P(M,s) 



4.An,b 



N-l i-oo 2 N 



X det[M,, jv(ejv)] det[Mi,jv(ejv)]5j(s)5i(s) , (CIS) 
where 



as given in the text in Eq. ( 95 ) . 



2. Non-intersecting Brownian excursions 

Li this subsection we compute the constant An,e of 
the formula (123) with the same method. Let us recall 
the joint distribution of M and tm for the N vicious 
excursions 



JV-l ( oo 



Pn^e{M,tm) 



5,(.) = lim/ d,?^^?fcile-^l^l 



]\,fN{2N+l)+3 



g2 + 2M2s 



(CM) 



'i— 1 n^ — 1 J 

X nN{{n,}\M,TM)nNi{n,}\M, 1 - tm), (C20) 
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and where 

oo 

f^Ar({n,}|M,0 = ^(-l)'=fc2A^(K},fc2)e 



k=l 



(C21) 

with the same notation as in the main text {n.i\ = 
(ni, . . . , nAf_i) and = (n^, . . . , n'j^_^). The normal- 
ization of this joint pdf is 

/•oo f 1 

/ dM dTM Pn,e{M,tm)^1- (C22) 
Jo Jo 

As in the case of bridges we can use the result of Ref. [50] 
of the cumulative distribution of the maximum of N vi- 
cious excursions 



FnMM) 



oo oo f N 



ni— 1 riN—l K. i—1 



(C23) 



where A?r is the square of the Vandermonde determinant 



AM{nl...,n%)^ det n'^^' '\ 



(C24) 



and ajf^E is a normalization constant obtained by impos- 
ing limA/^oo Fn,e{M) = 1: 



2iv^-jv/2j|^ir(2 + j)r(; 



(C25) 



The derivative of the cumulative gives the pdf of the max- 
imum 



PnMM) 



dFN,E{M) 



dM 



Til 



N 



(C26) 



The term ^p-riAr^ in the parenthesis of the second line 
comes from the derivative of the exponential, giving a 
factor ^ J2i '^i^ that can be simplified by the symmetry 
in the exchange of rii and njv- 

Now we have to compute the marginal Pn.e{M) by 
integrating the joint pdf given by Eq ( C20 ) over tm 



Pn.e{M)^ [ dTM Pn,e{M,tm) (C27) 

JQ 

to identify our constant An.e- As in the Bridges case, 



we compute instead the function of T 

N-l ( oo 



P(M, T) 



A 



N.E 



J^,.[N(2N+l)+-i 



X / dTM (K}|M, Tm) (K}|M, T - tm) , 
Jo 

(C28) 



and at the end we will recover 

PnMM)^P{M,T=1) 



(C29) 



The formula for P{M,T) is best treated with a Laplace 
transform, because the convolution in the second line of 



Eq. ( C28 ) becomes a simple product. Taking the Laplace 



transform with respect to T, one has 



P(M,s) 



dT e""^ P(M,r) 



A 



N-l ( oo 



N.E 



rife 2A/^ 



j^N{2N+l)+3 HE 

i=l (.ni = l } 

X (nNi{n,}\M,s)y , 
the Laplace transform of il^v being 

poo 

nN{{n,}\M, s)^ dT e-'^ 17jv({nJ|M, T) 
Jo 

^,fc2A7v(7^l^...,n^r_l2,fc2) 



(C30) 



= lim Y^i-e) 



>ik=i 



2M2 



fc2 



(C31) 



This e < 1 has been introduced to regularize the in- 
tegrand so that one can permute the integral from the 
Laplace transform and the infinite sum over k. Writting 
the Vandermonde determinant as 



N-l 



Ajy{ni^, . . . ,nAr_i^/c^) 



2 2\ 



X n (C32) 

l<i<j<A'-l 



one can use the identity, shown below in Eq. (Dl), with 
p2 ^ 2M'^s/tt'^ with s > to get 



nNi{n,}\M,s) = {-lYW2s 



M- 



nfj^ ( 



71-2 sinh {My/2s) 

X n ( 

l<i<j<JV-l 



(C33) 



Under this form, we are ready to take the inverse Laplace 
transform of P(M, s) using the Bromwich integral: 



P(Af,r) = 



1 

2i7r 



7+200 



ds e''^PiM,s), 



(C34) 



"f — tOO 
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where 7 is such that the vertical hne keeps all the poles 
of the integrand to its left. The poles in the complex 
plane are located in s — — (fc^7r^)/(2M^), for A; = 1, 2, . . . 
on the negative real axis, so that 7 = is well-suited. 
We close the contour with a semi-circle to the left. By 
the Jordan lemma, the integral on the semi-circle when 
i?, — > 00 is zero. The residue theorem then gives 



P(Af,T) 



A 



N-l ( 00 



N,E 



MN{2N+ 



k = l 



N-1 



l<i<j<N-l 1=1 



3 TT^ 

2 2A'P 



N-l 



2fc2 



(C35) 



Taking T = 1 and re-labelling k = n^, one obtains 



Pn.e{M)^P{M,T = 1 

00 00 N 



A 



N.E 



]l.fN{2N+l} + l \ ^2 



ni—1 njv— 1 i—1 



2 2M2 



N-l 

E 



2< 



(C36) 



In order to compare this formula with Eq. ( C26 1, we need 



to simplify the last term in the brackets. To this end, we 
introduce the function g{n) = g{ni, • • • , njv): 



N 



7in)^llnjA%{nl...,n%) 



(C37) 



which is totally symmetric under the exchange of any 
couple of its arguments. One has 



E--E E 



2n% 



9{n) 



1^1 f^N \ i—1 



E--E E 

rii riN \ i=l 



-1 "JV 



2« - nf) + 2nf 



N-l 



=E---E h(^-i) + E 



2n 



n 1 71 N 



N-l 



^E---E 2(^-i) + E 



2n 



gin), (C38) 



'AT 



where we change the role of rii and riAf in the last line. 



Equating the last line with the first, one obtains 



^N-l 



E-E E 



2n 



N 



5H = E---E(^-i)5H- 



9 9 

nij — ni , 

Til "N \ i=l JV I y „j 

(C39) 

Inserting this identity in formula (C36), the comparison 



with Eq. ( C26 1 is straightforward, and one obtains 



An,e iVvr^aAT.. 



2^^-iv/2j|V^-ir(2 + j)r(|+j) 



(C40) 



as given in the text in Eq. ( 120 ) 



Appendix D: Identity 



We want to show the identity 



sinh(7rp) 



(Dl) 



This will be shown by induction on N: we thus first 
analyse the case iV = 1 of this identity 



lim J^i-eY 



« ^-l fc=l 



p'^ + k^ 2 sinh(7rp) ' 



(D2) 



for any real p ^ (indeed this has a limit for p ~ and 
one can extend the result for p = 0). The first step is to 
separate the series in the Ihs into two parts 

00 00 9 

k P 



e — >lk=l e — >T-k=l ^ 



(D3) 



The first series is geometric and one has 

hm f;(-.)'^- = hm - 1) = -2 ■ (D4) 



e — >i k=l 

In the second, one can permute the limit with the sum 
because of the normal convergence 



lim 

e — ^1 fc=l 



p^ + k' 



00 

^ Z^^ ^ fc2+ 2 

fc=i ^ 

\k— — oo 

9 ' B + p2 2 



fc — — 00 



(D5) 
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To compute the rhs, we consider the function on the com- 
plex plane 



2 iiTZ 



1 



(D6) 



This function has simple poles in z = fc with fc G Z and 
in ±ip. One can verify that with z — i?e*^ 



1 



2|i?2 -p2| |sinh(7ri?sin6')| 



0. (D7) 



Then using the Jordan lemma, one has that the contour 
integral of f{z) along the infinite circle centered in z = 
is zero. Using the residue theorem, one finds 



feez 

+ {2in) 
+ {2in) 



2tT J fc2 + p2 

2ip e-2"P - 1 
-2ip e^^P - 1 ' 



or 



k p 



np 



fc2 + j)2 sinh(7rp) 



(D8) 



(D9) 



Inserting Eq. (|D9|) in Eq. ((DSj) and Eq. in Eq. (|D3|, 

one has 



+ 2 (,2sinh(7rp) 2 

2 sinh(7rp) 



1 / 1 7rp 



DIO) 



which is the announced result for iV = 1. 

Let us prove the general result by induction. Consider 
the quantity 



lim j2i-^y 

€ >1 k^l 

OO 

= lim ^(-e) 

e — >i k=l 



p2 -|- fc2 



and the induction hypothesis at rank N (Dl| that we 



use in the second term of Eq. (Dll I, we have proved the 



formula at rank A'^ + 1 which ends the induction proof of 
the identity ( [Pl] ). 

In order to prove the lemma (D12|, we put e = e^^ 



with the real x > 0, and we shall take the limit x 
Denoting for q = 0,1,2, .. . 



0. 



one sees that 



k = l 



~^9qix) = gq+l{x). 



(D13) 



(D14) 



Using the fact that go{x) is a simple geometric series, we 
apply this relation to compute gi{x): 



dx 

i! 

da;2 Vl 



e-^ (e-^ - 1) 



(1 + e-^)' 



(D15) 



The domain of definition of gi{x) is the positive real axis, 
and in this domain, 51(2^) coincides with the function 
defined over all the real line 



X gi{x) = 



e ^ (e 



1) 



(D16) 



p2 + fc2 



JV-1 



= lim^(-ef e\{{e-n^) 



One sees easily that gi{x) is an odd function gi{—x) ~ 
—gi{x), which implies that limj-^o 51(2;) = 5i(0) 0, 
which shows the lemma for q — 1. To go further, the 
Taylor series of gi{x) must have only odd powers of x. 
Because they coincide on a; > 0, this is true for 171(2;): 

00 

V2;>0, .gi(x) ^ai,pa;2p+i. (D17) 

J9=0 

Indeed this odd power series expansion extends to 
all g q{x) (with other coefficients) thanks to the rela- 
tion (D14) because we differentiate twice to obtain gq+i 
from gq 



e — ^1 k=l 



lim(p^-Kn.^)f:(-.)^ ^^nfy;(fc^-n.-) 



i=l 



Vx > 0, gq[x) = Q(j,p2:^^+\ 

p=0 



k=l 



With the lemma (which we show below) 
00 

lim ^(-e)'= /fc29 ^ 0, VgeN*, 



(Dll) 



(D12) 



and this prove that 



lim gq{x) = 0, 



(D18) 



(D19) 



c — >'l k=l 



which, regarding the definition (D13), ends the proof of 
the lemma (D12 1. 
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Appendix E: Some useful formulae for Hermite 
polynomials 



For the stars configurations we use tlie identity 



For our purpose, the following integral representation 
of the Hermite polynomial Hn{z) is useful 



2"+! ^ foo 



dte~* cos {2zt-n7r/2) (El) 



The formula used in the text is explicitely 



/■oo ^ 

J dq q'^ cos [qC - n- 



M 



n+l 



M 



(E3) 



dq g" cos ( q — n— ^ e 2m^'^ 



2t 



— ili 7-7- 

e 2t iy„ 



7\f 
2i 



(E2) 



(in particular with t = I — tm and n — 1). 
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